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Steady State Solution of the Relativistic Boltzmann 
Transport Equation’) 


Ivan ABponyl, Budapest, Hungary?) 


Introduction 


A detailed study of the Maxwellian velocity distribution function was made 
by F. JUTTNER [1]’) as early as in 1911 from the point of view of relativistic 
mechanics. JUTTNER’S treatment is based on a generalisation of the familiar 
phase-space method in which the Liouville theorem holds for the motions 
described by the relativistic equations of motion [2]. Then the distribution 
function is determined in the usual combinatorial way. 

Here we propose to follow another line of thought. Namely, we start from the 
relativistic Boltzmann transport equation given by CLEMMow and WILLSON [3] 
and KeEn-1TI Goto [4] and complete it with the collision integral. Then we try 
to find the steady state solution of this equation. 


The Relativistic Boltzmann Transport Equation 


Let f(t, r,u) denote the distribution function of a system consisting of N 
identical particles, r the position and 


fe nee ETS (1) 
v \2 
n- © 
the reduced velocity of the particles (c being the velocity of light in vacuo). 
We note that 


v \2)-1/2 i 9) 1/2 
ati lata eee): 
The particles of rest mass m, obey the relativistic equations of motion 
7 aa F ? (2) 


} 1) A part of these calculations was carried out while the author was the guest of the Osservatorio 
- Astrofisico di Arcetri, Firenze whose kind hospitality is gratefully acknowledged. 
¢ 2) Institute for Theoretical Physics, E6tvés University. 

3) Numbers in brackets refer to References, page 175. 
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where F is the ordinary force vector. The relativistic Boltzmann transport 
equation as given by CLEMMON and WILLSON is 


At 
af tas VII ae EV + et | 2s) 


where ,, is the gradient operating on the functions of u, (Af/A?),,,;, represents 
the change in the value of # produced by collisions during a time interval of 
unit length (invariant collision integral). CLEmmMow and WILLSON neglect the 
expression (Af/At),,1;, because they do not take into account the role of the 
collisions. Though their approximation is appropriate for many practical 
applications, a more precise treatment of certain theoretical and practical 
problems demands closer consideration of this term. 

As is well known from text-books of statistical mechanics [5], for the mo- 
tions described by the Newtonian equations of motion the collision integral 
can be written in the form 


ae A =|[rr ‘) g L(g, 8) dQ dv’, (4) 


where } = f(é, r, v) is the distribution function for the system, and 


f=itrne), f=ftro, f=tero), (5) 


v and v’, v and v’ being the respective velocities of the colliding particles long 
before and long after the encounter; I(g, #) is the differential cross section of 
the collision; and 

g=|v—v'|. 


Remembering the deduction of (4), it is easy to see that (4) is valid also in 
the case of relativistic motions described by (2). The only change to be made 
is the substitution for the ordinary velocity v of the reduced velocity u. Then 
the relativistic form of the collision integral is 


(Sr, = [CF - #1) ete, 8) dQ aw’, (6) 


where g = |u—u’'|, and 


f=tire), f=ttrne), fier, f=fenw), My 


u and uw’, u and uw’ being the respective reduced velocities of the particles long 


before iil long after the collisions. 
CLEmMow and WILLson have pointed out that the left hand side of (3) is 


an invariant scalar with respect to the Lorentz transformations. It is easy to 


see also that (6) is an invariant scalar. 
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So the complete Boltzmann transport equation in relativistic form looks 
like 


of G 
oft 


vital (uV) f+ me {(FV,) + (V, F)}/ 
ae (8) 
=| (¢f —ff) g I(g, 9) dQ du’ . 


In general — even in the case of the velocity-dependent Lorentz force — 
the term containing 


VF 


makes no contribution, and we shall therefore omit it in the following. 


Determination of the Steady State Solution 
We may call a state steady when 


of 
sae 


and when the collision integral has no effect on the distribution; in other words, 
when the state is not changed either by the time or by the collisions. These two 
suppositions made, we have to solve the differential equation 


c 
yit+w 


Later we shall have to demonstrate, that the solution of (9) is consistent with 
the second supposition, i.e., that the solution makes the collision integral zero. 

For the time being we shall confine ourselves to one spatial direction and 
conservative velocity-independent forces only. Then the differential equation 
to be solved is 


(UV) f+ a> (FV) f= 9. (9) 


cu (L +m)? 9 + (ong 0) F(x) 2 = 0. (10) 


This can be separated. Supposing a solution of the form 


f(x, u) = A(x) > Blu) 
we get 


Date ah bg th tt Bily) 


F(x) ~ A(x) um,c? Bu) nde 
from which we easily compute that 


a, a= oO exp{— aM, Cc? y1 + u*—a V(x)} : 
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Here V(x) == f F( F(x) dx, and a and 0 are the constants of integration to be 
er later. In ie general case we find that 


jr, u)=6 exp{— aM, Cc? yi+u +u?—a V(n)} : (11) 


We can compute S, the entropy of the system under consideration, with the 
help of the distribution function in the form 


Core h ff te. u) log f(r, u) dr du, 


where k is the Boltzmann constant and the integration is extended over the 
whole phase space. Then denoting 


N = | He. u) dr du, (12) 


E=N B= | {me ct + u? + V(n)} f(r, u) dr du, (13) 
we obtain 
S=—kNlogb+akE. 


Taking into account the well-known thermodynamic ‘relationship between 
energy EF, entropy S and absolute temperature T: 


Os ee! 
(oe), = 7 


1 
a= a. (14) 


we can conclude that 


In order to see the meaning of the other constant of integration we consider 
that we are describing a system consisting of N identical particles. The total 
number of particles, in the case of zero potential energy, is given by 


N= [/ ir, dr du = bY [J exo| 96 Tee du, 


where V is the volume of the system. Performing a polar transformation of the 
variables of integration, we arrive at 


IN ‘aa fexol 2 [1+ ah ut de, 
u being the length of the vector u. Let 


; 2 
a= eT. and “= sh@ ‘ 


——— 
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then 


N=426b v [esol a chp} sh? ch@ dp 
= 2bvV fest a che} (ch 3 y — ch) dp. 


This integral can be evaluated with the help of the following formula [6]: 


fo,e) 


[expt z cht} chp ¢ dt = —— e-? (4?) F(z) (15a) 


0 


which is valid for 0 S argz <a and p > 0. Here H!) is the Hankel function 
of the first kind with the index #, obeying the following recurrence formulae: 


2 Hp” (2) + 2 Ap (2) = 2p Hp (2), 


(1) (1) dy oo) 
Hy” 1 (2) — Ay (2 y= 2 vy HT, (z) . 
So N can be computed to be 
N= dV (HM) + H(@)}= 22°" HQ). (16) 
Then 
j Ce 
N a N ah 
b= — (22?)- —— (22?) ‘ 17 
A Goins tae semen re reap oe 


And so the distribution function of the steady state is 


5 oe 
N kT 


Oem aeeraces 


exp| -# (mg ee y1 + 2 + vin))}. (18) 


It is easy to see that the collision integral (6) vanishes for this solution. 
One can demonstrate this statement by direct substitution, taking into account 
the fact that in the collision the normal components of the reduced velocities 
_ change their sign. So we can conclude that this solution is the appropriate one. 

In general, we have to use the complicated expression (18) which can be 
evaluated for a given value of m, c?/kT by using tables of the Hankel functions. 
But in the limiting case, when m,c?/kT >1 we can obtain explicit results 
without such difficulties, making use of the asymptotic form of the Hankel 
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functions valid for large values of the argument: 


: . Z, 
HO (i x) = 1 Var ae (19) 


x being real and positive. Then we obtain instead of (17) 


N Mo C* \3/2 My Cc? 
oe owe! exP| er} (20) 
and so the distribution function in this approximation becomes 
N My C* \3/2 My Cc ones 
i(r,u) = + (seer) exp|— Fe (ee 1}. (21) 


This may be called the relativistic Maxwellian velocity distribution function. 
In order to see its connection with the non-relativistic function, we transform 
the expression 
f(r, u) dr du 


into a function of r and v. Evidently 


j(r,v) dr do = f(r, u) ota dndes 


where f(r, v) is the distribution function with the new arguments and 0(u)/0(v) 
the Jacobian of the substitution. Taking into account that 


So we see that 


fir) dr do = F(z)" [1 > (2) 9] exp 


which is exactly that given by JiTTNeER [1]. 


Remarks 


_ Here we may call the attention of the reader to an already well-known 
result, namely, that in general the law of the equipartition of the energy loses 
its meaning in relativistic statistics, because we cannot write the a 


co ae 
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kinetic-energy as a homogeneous quadratic fotm of the velocity variables. In 
an unspecified case of mg c?/k T with zero potential energy, using (13), (15a) and 
(15b), one obtains the average total energy of the particle in the form 


‘ . Hy’ (i ( “a 
E=kT{1 bi Lat ap ( mm) 


(dash denotes differentiation with respect to the whole argument). Even in the 
case My C2 > k T this expression gives 


Ba kT SMe. 
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Résumé 


Depuis les travaux de JUTTNER [1], nous connaissons la forme relativiste de la 
fonction de répartition maxwellienne des vitesses. Le traité de JUTTNER est basé 
sur une généralisation de la méthode de l’extension en phases, de sorte que le 
théoréme de Liouville est valable par rapport aux équations relativistes de 
mouvement. La fonction de répartition est déterminée par la méthode combinatoire 
usuelle. 

Ici, nous proposons de suivre une autre voie. Partant d’une forme relativiste de 
l’équation de transport de Boltzmann donnée par CLEmMow et WILLSON [3] et 
Ken-i111 Goro [4], nous proposons une expression relativiste de l’intégrale de col- 
lision qui compléte cette équation, et nous donnons ensuite sa solution pour l'état 
stationnaire. 


(Received: August 21, 1959.) 
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The Transverse Flexure of a Thin Circular Plate 
Subject to Parabolic Loading over a Concentric Ellipse 
By W. A. Bassatt, Alexandria, Egypt’) 


1. Introduction 


The behaviour of transversely loaded elastic plates has for several years been 
the subject of much research. Many analyses have appeared of cases where the 
boundary of the plate is rigidly clamped, simply supported, or free, and the 
plate is under various distributions of normal thrust; excellent bibliographies 
are to be found in MUSKHELISHVILI’s book [1]#) and at the end of a recent 
paper by DEveERALL [2]. The problem has been extensively investigated in a 
series of papers [3, 4, 5, 6] for a thin circular plate subject to different types of 
loading over a concentric or an eccentric circle, when the edge of the plate is 


elastically restrained according to a general boundary condition which includes ~ 


the rigidly clamped and simply supported boundaries as particular cases. When 
the boundary of the circular plate is free, and the plate is subject to sym- 
metrical or uniformly varying pressures over an eccentric circle and is supported 
at several interior or boundary points, exact solutions for the problem have 
been developed by the author [7, 8] through the application of MUSKHELI- 
SHVILI’S complex variable method. Appropriate expressions for the complex 
potentials, deflections and stresses have also been established [9, 10] for an 
elastically restrained circular plate which is partially loaded over a sector or 
over a segment of the plate. In two recent papers [11, 12] the bending of the 
same plate has been discussed, when uniform or hydrostatic loadings are 
distributed over the surface of an ellipse concentric with the plate. Making use 
of MUSKHELISHVILI’s complex variable solution of the first boundary-value 
problem in the theory of elasticity for elliptic plates ([1], p. 244), SaPONDZYAN 
[13] has determined the transverse deflection of a clamped elliptic plate uni- 
formly loaded over a concentric circle or over a confocal ellipse. 

In this paper we follow the same procedures adopted in the two previous 
papers [11,12], in order to derive the complex potentials which solve the 
problem of an elastically restrained circular plate subject to parabolic loading 
extending over the area of a concentric ellipse. Standard formulae for the 


complex combinations of stress components are used in the computation of the — 


*) Department of Applied Mathematics, Faculty of Science, University of Alexandria. 
*) Numbers in brackets refer to References, page 190. 


——— ee ee 
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boundary and central values of the bending, twisting moments and shearing 
forces. The two limiting cases in which the loaded elliptic patch becomes a 
circle or a line ellipse are also considered. 


2. Basic Equations 


Let C denote the boundary of a thin circular plate, of centre O and radius c, 
made of isotropic homogeneous material whose Potsson’s ratio is 1 and 
Youne’s modulus is E. If 2h is the thickness of the plate then its flexural 
rigidity is known to be given ((14}, p. 3) by 

h> 
Da2 E 30-7)" (2.1) 
Let z=x+iy=re'® be any point in the mid-plane of the plate which is 
chosen as the plane Z = 0 of a rectangular Cartesian frame O(x, y, z). We 
assume that the plate is divided into the two regions 1 and 2 by the contour I" 
of the ellipse 


é (2. 2a) 


and we write 
@=2+?, P=@-—-P=aae, (2. 2b) 


where ¢ is the eccentricity of the ellipse (see figure). If £,(z, z) (j = 1, 2) are the 


transverse load intensities at points of regions 1 (inside J’) and 2 (between I” 
and C) then, according to the Poisson-Kirchhoff theory of the bending of thin 
plates ([14], p. 88), the small transverse displacements w, ( =1, 2), measured 
in the upward direction OZ, at the point z satisfy the Lagrange partial diffe- 
rential equations?) 

DV4 w, (z, 2) = — (2, 2) G= 1,2), (2.3) 


3) Bars are used to denote conjugate complex quantities. 
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where 2 denotes the Laplacian operator 02/0x* + 07/0y? = 4 0?/ 0z 0z. The 
general solutions of equations (2.3) may be expressed in terms of four analytic 
functions Q,(z), «;(z) (j = 1, 2), usually called the complex potentials and two 
particular integrals W,(z, z) corresponding to the load intensities ;(z, 2): 


w, = 29,2) +232) +o) +5 +W22) G=12). 4) 


Adopting STEVENSON’s nomenclature*) and method of specification of 
stress [15] the bending, twisting moments and shearing forces at any point z 
of the plate are furnished by 


10 — Or = —4 (1 +n) P [2 Re Q(e) + 552 


Oz 0Z Aisa 


~ ~~ ~ 2 wv Ww; 

76 + 6r — 2160 = —4(1—7) P [22% Bs ESC aay (2. 5b) 
~~ — 3TV. 
(Zi 62 = —8 P(*) [O") + ae (2. 5c) 


where the plate modulus P denotes Love’s flexural rigidity D divided by the 
thickness of the plate, Re stands for the real part and the primes on Q,, w,, 
and W, designate differentiation with respect to z. 

The transition conditions expressing the continuity of the deflections, 
slopes, moments and shears at any point on I’ may be taken as 


feb [eel eae 20 


These conditions have been found extremely helpful in treating problems 


concerning partially and normally loaded thin circular plates (see References 
[3-12]). 


3. Boundary Conditions and Boundary Values of Moments and Shears 


Along the circular boundary C of the plate we shall assume that 


v6 +4 
v= 0 ) ao ¥ ’ (3.1) 
Ov B-—A 
4) The relations between STEVENSON’s notations [15] and those used by TimosHEnko [14] are 
M,=—2hns, Ms=2hsn, Mys= —Mgn =2h5S = — 2htn, 
Qn =- QhnZ, Q,= —2hsZ, 


2h being the thickness of the plate. 
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where A is a dimensionless restraining parameter and f = (1 — y)/(1 +). It 
was proved [3, 9] that the conditions (3.1) imply that 


0?w Ow 
w=0, € ot Bag eS = 7, (3:2) 
OPW, 
n—? yvy=1° 6 Oy a as. ee) 
60 SbF a ak Ow 

ten 4 y 6. Or 0e-* (3.4) 

,. OW, (or, Slob t+ WD) 
Re [225 + 552 — A282, ee © “|| =0 (3.5) 

along C, where 

oe! Step ea 1 Si 3.6 
ai ea eRe we la Ae —p (3.6) 


The boundary C is clamped or simply supported according as A = 1 (vy = oo) 
or A= —f (y=™%) , respectively. Equations (3.3) and (3.4) are convenient in 
computing the boundary values of the moments and shearing force 6Z. The 


boundary value of the shearing force rZ can be calculated by using (2.5c). All 
these equations have the advantage of not involving the complex potential 
(2). 


4. Complex Potentials for a Thin Circular Plate under Parabolic 
Loading over a Concentric Ellipse 


In two previous papers [11,12] the particular integrals and complex 
potentials corresponding to the load intensities p,(z, z) = constant po, p, = 0 
and ,(z, 2) = py X = by (2 + 2)/2, Pp =O have been completely determined. 
We now deal with the case 


Pilz, 2) = Po = fo 22, Polz,2) = 0. (4.1) 


In this case the total load on the plate is easily seen to be x py a b d?/4 and the 
particular integrals may be taken as 


RO? W,{2,2) =O, (4.2) 
where 
he Po ee Po (4. 3) 
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It is known that the conformal transformation 


z= FZ [lato o+ (a— 07] (4.4) 


maps J’ on the boundary of the unit circle in the ¢-plane. The point ¢ corre- 
sponding to the point z on J’is therefore given by 
ZZ a B= 


a Pe (4.5) 


eae a—b’ 


where?) 
age eee CP (4.6) 


Substituting from (4.2) in (2.4) and introducing the resulting expressions in 
the continuity conditions (2.6) we find, using (4.4) and noticing that along J’ 


we have fe 25, (G3 


[22(2) + 22@ +02) +O@k= 7 kA P= eke +P (+ CP, 


(4.7) 

[z Q'(z) + Q(z) + w'(z)] | 
it = (4.8) 
=ZhePa sekla-HNC+@+HNCN RHP +Cy, | > 
(Qe) +O@k=kAP=T ke EPH EP, 4.9) 


[Q"() = Zhe = h[(a—d) C+ (a +8) C4] (2d? + (2+ C-%)]. (4.10) 


From (4.4) we deduce that 


(4.11) 
and hence (4.10) gives 
/ il = 
[ae 2 O) =F bab+ PC -—Cy Ree P+ ey. (4.12) 
Integrating (4.12) and introducing for convenience the notations 


Up te Co yea) ee eaten (4.13) 
we immediately get 


Oi) eae Flea bal ln ee 
Moly. ob 4 ffm t Pd uy + 2abf+ Ky], (4.14) 


°) No confusion need arise between this definition of Z and the third dimension of the pit 


as expressed in the stresses Z, 
ag ih 02, XZ, and yZ, since the third dimension Z only arises in such 
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a 


where K, is a real constant*) to be determined. Substituting from (4.14) in 
(4.9) we conclude that K, = d? (3 a® + b?)/2, and hence, using (4.11), equation 
(4.14) becomes 


[Q"(z)]3 - [= ae ee a =%) é-2 a 2(2)| 


1 = 1 
=k [sa bd® loge + 7 ftuy + fd? u,+ 2ab fr, (4.15) 
tl 
+5 @ (3 at+ 0) 
The hee of (4.15) yields 
(2 (k= 5 | g b [16.4 bd? zloge +5 ae (a uz + b v,) 
+ oe f? {a (3 a® + 13 6?) ug + 6b (13 a? + b?) v4} 


, (4.16) 
+ 5 4 (at — 8 a? — 17 4) u 
— 5 b (15 at + 8 a2 0? + bf) vy + Kal, 


where K, is an arbitrary constant that may be complex. From (4.4), (4.15) and 
(4.16) it is found that 


[2.0'(2) + O@E= ze k [ay f(a ms — 6 v5) 
+ = PP {a (9 a? — BY) uy + b (a2 — 90%) vg} ¢ (4.17) 
+a (3at— 1104) m+ b Vat —b4) v, + Ki. 


Subtracting (4.17) from (4.7) leads to 


, L 2 dad 1 
fo'(z)le = & Peele lhc? ae o(0)], 
jel FG Bo a i ans 
=, kf {3 (bv; — @Us) — a Us 3 (4.18) 
+ a (17 bf + a? b? — 2 at) wy ; 
— b (b4 + a? b? + 14 at) v K,],| 
8) A pure imaginary constant 7K in [Q’(z)4 contributes a term iK z to Q,(z) or Q,(z) and can 


therefore be dropped since it makes no contribution to w. 
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which, on integration, gives 
[oleh = ay [4b (6at + a? b? + BA) log’ 
— {5 (at + 17 a® BD? — 2 4) uy + Bab d? 0, (4.19) 


1 1 
— =f (@u+3abu%)—z PF 


where K, is a constant of integration. Since a pure imaginary constant in «(z) 
leads to w = 0 everywhere it follows that K, may be taken as real. Inserting 
(4.16) and (4.19) in the remaining transition condition (4.7) we deduce that 


Ky = % @ (7 ¢ + 107 a? b — 2 as). (4.20) 


From (2.4) it follows that the terms 


i 


6 * Ke and ah Kya 


entering in (4.16) and (4.19) contribute nothing to w and can accordingly be 
ignored. Introducing the value of K, from (4.20) in (4.19) we finally see that 
the transition conditions along J" are satisfied by taking 


1 
[Q()B = =p h [16 abd? zlogl — 5 a (17 4 +8 a? B* — at) wy 


1 
— = 5 (15 at + 8 a? b2 + OA) v, 


1 (4.21) 
qx {aBa+13 b?) us + b (13 a2 + 3 B2) vg 
1 

+ 20 (a us + bv5)], 

[w(z)], = gg b [4b Oat + a2 5 +9 04) loge 

1 

+] 4 (7 b4 + 107 a? b2 — 2 at) 

(4.22) 


—fi> (a4 +17 ab? — 208) w+ 8abd* v4} 


ed 


-~Fh@ut3abv)—2 9 ug] 


ote Ome « 
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Substituting for ¢ and ¢-! from (4.5) in (4.21), (4.22) and remembering the 
definitions of u, and v, given in (4.13) we obtain, after considerable reduction 


ee z+Z 
[Q2)k= hab [Fe +e (2 log= 5 — Z Gz |, (4.23) 
[w(2)]h = kab [H(2) + 4 a® 7? {(16 « — 19 +9 73) jog te Z K(c 2), 
(4.24) 
where 
Z= (at — faye, <r = 6-2, (4.25) 
4 22 4 
Fee) =ab2 [1+ ;at-3% + (Ba—-8 +0 vara | a Ze 
il uh 6 24 
G(2) = +5 (16 = oak (4.27) 
HQ) = 70b@ [1-7 + 3-4 (6a 16 +a) BS, (4.28) 
K(2) = (s52) [48% — 63 + 49a + 4 (88.4 — 88+ 9 or) i 
a (4.29) 
— 16 (16 « — 16 + 3.0) ule | 


Letting b tend to a and « tend to zero in equations (4.21), (4.22) or equations 
(4.23), (4.24) we find that they simplify considerably to 


[Oleh => kz at [4 log= — 3], (4. 30) 
lo(2)]i = — #a® [6 log +7]. (4.31) 


It is now immediately checked that (4.30) and (4.31) are in agreement with 
equations (2.13) and (2.16) of [7], p. 732 in the special case % =0, n= 4. 
J We now introduce the notations’) 


C2 


ae Z' = (f'? — 22)U2, (4.32) 


7) Accents in equation (4.32) do not indicate differentiation with respect to 2. 
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Assuming that ” and s are any two integers including zero it can be easily 
proved that along C we have 
i n Ze 5 
rel) GI: 


Re [(7)' (7)] 

Re [2 Z (3) = Re [z 2’ (Ey: (4.33) 

ee) 
zZ+Z 


I 


Since the expressions involving Z in equations (4.23), (4.24) are not uniform’) 
inside I” (region 1) we deduce that these expressions should enter into Q,(z), 
@,(z), respectively and not in 2,(z), w,(z). Moreover, since we require w, to 
vanish on C it follows that these expressions containing Z have to be associated 
with other terms involving Z’ such that Re [z 2,(z) + w,(z)] vanishes along C. 
Equations (4.33) are applied in the derivation of the appropriate expressions to 
be added and the continuity conditions imply that they should appear in 
22, (2), @,(z) as well as in Q,(z), w.(z). It is found that these additional expressions 
containing Z’ have poles at z = 0. Terms involving negative powers of z have 
therefore to be added in 2,(z), w,(z) in order to remove these singularities at 
z= 0. To satisfy the transition conditions, these negative power terms must 
also occur in {2,(z) and @,(z). Equations (4.33) are then used to combine them 
with the appropriate positive powers of z so as to get zero deflection on C. 
Carrying out the foregoing procedure we find that the appropriate forms for 
the complex potentials 2,(z), w,(z) may now be written as 


O,2) = kad [x g(2) + a2 fz jog 50 4 ZG(z) — 2’ L(2) 


PI 4%, (4.34) 
Ts i i) 3 (5 a aie 
@(z) = kab E c® p(z) + iz a? P {(16 a— 19+ 94-4) log= 7 | 
+ ZK(2) — Z' T(2) +3 (304 — 1) . = =) 
a i / (4.35) 
+4(8a—-84a) (— 4) | 
— 2 (16.0 — 16 +347) (S - iI, 


*) The two branches of Z = (2? — f2) 


12. interchange wh i ; . 
focus of the ellipse I’. g en 2 describes a closed path round either 


om Bowe 
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where ¢(z) is a function regular within region 2 and 


wi tiga ae Gps? 
FAS) rs ( ee ee ), see) 
T(z) = a [48 a — 63 + 49 a4 + 4 (88% — 88 + 94-1) - 
(4.37) 


ve 
es 

. 
— 


~ 16 (16 a—16 + 3a) Fo 
Adding (4.34) to (4.23) and (4.35) to (4.24) we get 


Q,(2) = kad [z (2) + F(z) + as Dees, et — Z' L(2) 


a hicis ons (4.38) 
eles eaeneensy ye 
e(e) = had [— c gle) + He) + gy a? {(16 0 — 19 + 94) log FEF 
— Z'T (2) +3(3a2—1) (G = =) 44 (8a—8 +07) (4.39) 
x (& be ‘a ) — <2 (16 « — 16 + 3.0) (& - 4)\]. 


Since the function y(z) appears in 2,(z), @,(z) as well as in 2,(z), w.(z) it follows 
that this function must be regular all over the entire plate. We have now to 
determine g(z) from the boundary condition (3.5) in which we put W, = 0. 
The following terminology has been introduced in a previous paper [12] and 
will be used here. If f(¢) is a meromorphic function of the complex variable 
t (|¢] <u <1) which has a single pole of order s at ¢ = 0 with principal part 


yes 
n=1 he 
then 
Reg f(t) = f() — gf) - (4.40) 
Adopting this notation and applying (4.33) it can be shown from equations 
(4.34)-(4.37) that along C we obtain after some algebraic manipulation 


Re Q(z) = kab Re [{zg'@)y’ + a2 {1 + 42 Reg og (2+ 5}, 4.40) 


Re [2Q"(2)] = kab Re [2 o'@} + Pi +5 + 42 Reg _ (2+ ae )}] , (4.42) 
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ie u I a f Oy ey 
Re[2 of(2)] =—kad Re [22 9"@) + ty 08 GG (160-19 + 9art} 
+ 4+ Regl% (640-79 + 47 | 
ip L (4.43) 
416 (124 — 12 + ot) Fe | 
Bite ie 
16 (160-16 +30) ZI]. | 
Inserting these expressions in (3.5) with W, = 0 leads to 
! 1 Z’3 3 f2 
Re[p + (1-4) z9'@) + (1-4) ae — > Rega (2+-5 )} 
A a® ie ; % ne Z 
24 c2 {16 a 19+ 9at + 2 = Reg (640 — 79 +4701 S (4.44) 


ie 
BStO (120 12" at) 10 oie = 16 ees) HI sis 


on C. This can be satisfied by equating the expression between the square 


brackets to i. K, where K, is a real constant that can be ignored since it makes — 


no contribution to w. Introducing for convenience 
’ 
De p(z) = pit), de me aay en oe.) (4.45) 


we get the linear differential equation 


: 
ty'(t) + my(t)=y+ De y, Reg iy & (l= 72 )ua)s (4.46) © 
7=0 
where 
iA G41) de (16 —19¢@14 9q-2 F 
2 24 a +Ie ere | 
16 
ore (m + 1) d? 


1 1 

Vi = se [8 (mm + 1) d2 — 5 (m= 1) (64 — 79.02 + 47 0-2) at 47) 
8 

yam —5 [(m +1) d+ (m1) (4$—4et4 4 a*) SI, 


8 4 
Ve= qq (m—1) (16-1604 + 30-2) 2, 


~~ 
a 
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For a clamped boundary 4 = 1, m= oo and the function ¢(z) is immediately 
obtained by putting A= 1 in the expression between the square brackets in 
(4.44) and equating the result to zero. 


In the general case where m + oo the solution of (4.46) is 


3 
vt) = 5 + Dy Si) + Kye, 


j-0 


where 


t 
Si(t) =t-™ / —l Reg [t-2/ (1 — 2)¥2) dt (7 = 0,1,2,3), ~ (4.48) 
0 


and K; is a constant of integration which, since m > 0, must be zero, so that 
w(t) becomes regular at #= 0. From (4.40) we have 


j-1 co 
Reg [¢-24 (1 — 22)4/2] = ¢-2i t — fie — dc, es = = Cn4," (720), (4.49) 
n=0 n= 


where 


@=1, 4=-=5, GE a ee = (n> 1). (4.50) 


Introducing (4.49) in (4.48) we obtain 


t 


Si(t) = ee pie [ (1 — #)22 — : oa el Git OLR Baron 
n=0 


0 
E where 
&= zt (j= 0, 1,2, 3). (4.52) 
We therefore have 
3 love) 
pt) = vt) =F + 4 0 = ay FD 8, 2", (4.53) 
7=0 n= 
where 
3 
Yo On + V1 n41 + V2 nt2 + Vs nts 54 
= ds 1 on = : Ton ¢ (4.54) 


The integral in (4.51) can be written in the equivalent form 


Si(t) =" fe cise at ot icy (1 —vjl do, (4.55) 


(1—#)1 
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and for even positive integral values of m it can be easily evaluated in terms of 


elementary functions. This corresponds to values of the restraining parameter re 


of the form 1 — n-! where 1 is an even positive integer. Having obtained the 
function g(z), the particular integrals (4.2) and the complex potentials (4.34)— 
(4.39) corresponding to the boundary conditions (3.2) and to the load (4.1) 
over the concentric ellipse are now completely determined. With these values 
inserted in (2.4) the deflection at any point of the plate may be computed. 


5. Central Values of Deflection, Moments and Shears 


Letting z tend to zero in (4.39) and (4.37) we find that the deflection at the 
centre of the plate is furnished by 


14(0) =2.0,(0) = 2k ab [H(0) — 2 g(0) + oN Se ie OBES 


: (5.1) 
+ (16 — 19 e? + 9 ¢4) logs || : 


Substituting for H(0) from (4.28) and for p(0) from (4.53), (4.54), (4.50), (4.47) 
and using the value of m given in (4.45) we obtain 


(ay Apt 20 = 
w,(0) = 25> [5 ats — ee + Set t (16-198 +94 


(5.2). 


x (log 7" 7} wes emceen 


Setting z = 0 in equations (2.5) leads to 


iiiy = Wy = 1Ly = yZy = 0, 
4Vq — V%q = —8 (1+ 7) PQ{(0), (5.3) 
#Yo + Pky = —4 (1 — n) Paw; (0). 


From (4.38) and (4. 3) it can be shown after some straightforward calculations 


that 
Q(0) = kad |a (1 1 1 
og ay 5A 5 +ab(1+— 2) 
A 4 
+, 3g (10 10 iens 
= 3 b d? 2 (A — 2) d® + (2 — 19 ¢2/8 4) g4/r2 
ey (0) = ha b(t ae ee ) ant seer. 


( 


With these values introduced in (5.3) the bending moments xy, and YXq are 


determined. 


“ay «* 
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6. Boundary Values of Slope, Moments and Shears 


Substituting from (4.41) in (3.3) and remembering that W,=0 and 
k = p,/(128 h P) we find the following formulae for the boundary values of the 
bending moments and slope: 


£ oo io ab 
Or gac AG he 1A 


aw net. 


v0 or We (2 
c 


7 (6.1) 
x Rel {z ple hd? + Fo dt Reg Fe (2+=5)], | 


where z = c e*®. Inserting the value of p(z) given by (4.53) we see that the right 
side of (6.1) may be written as 


Po ab 7 
CE one jer + + By 6, u2" cos2 n |, (6.2) 
where u = f/c, 
BS (2 $1) 8.4 Fe PB extn — usa = 200 (6.3) 


and y,c,,0, are given by (4.47), (4.50), (4.54), respectively. Differentiating 
* (6.2) with respect to @ gives 
ea 
c \ Or 06 Jy=c 


and hence from (3.4) we get for the boundary values of the twisting moment 60 


and the shearing force 0Z 


RA ee Se ee ae 
mate ayy ears ee fa eae sin2n@. (6.4) 


There remains the determination of the boundary value of the shearing force 
vZ by substituting from (4.42) in (2.5c) which yields 


; 16 z’3 3 f”? 
aay fa bar? Re | fz? ¢'( MY + a+ Ge @ Res (2 +75 )], (6.5) 
_ and substituting from (4.53) we obtain 
A ee ee ee G + 7 6! w2" cos2.n i), (6.6) 
n=0 
_ where 
8! =2n(2n 41) 6,4 ap PB 42 — San — 2%) (6.7) 


190 W. A. BassALt ZAMP 


7. Limiting Cases 


As a check on the validity of the foregoing results we deduce as a limiting 
case the complex potentials corresponding to the same boundary conditions 
(3.2) and the same loading (4.1) distributed over the area of a concentric 
circle. Letting 6 tend to a and « (= e~?) tend to infinity in equations (4.34)- 
(4.39) leads after reduction to 


ON ae E +A/ oo 5) log-“], (7.1) 

ox(e) = — 0% [5 4+ (1-54) +e loge|, 7.2) 

Oylz) = Po [1 +a( A ~ 5) — ee], (7.3) 
tog(e)ie ee [;4+ foe 2) ee log], (7.4) 


which agree with the expressions obtained by putting s = 4 and s’=6 in 
equations (3.39) and (3.40) of [3], page 15. Expressions for the central values — 
of the deflection and moments and the boundary values of the slope, moments | 
and shears immediately follow by letting 6 tend to a (e > 0) in the results of 
sections 5 and 6 and it is easily verified that the formulae obtained coincide 
with those given®) on page 597 of [4] in the special case n = 2. 

In the limiting case when the minor axis of the loaded elliptic patch tends 
to zero we have a line loading with variable load intensity extending along 
the real axis from x = — a to x =a. Assuming that b tends to zero and 9 
tends to co such that 2b f, tends to p, it is seen that the intensity # of this 
line loading at any point distant x from the centre is given by 

p=hi ge (1 a =)". 


az 


(7.5) 


The complex potentials appropriate to this line loading are immediately 
obtained by putting p,/2 instead of b 49, f =a, f’ = c/a, d® = a? anda =1in 
equations (4.34), (4.35), (4.27), (4.29), (4.53), (4.54), and (4.47). 
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Zusammenfassung 


Auf der Basis der Poisson-Kirehhoffschen Theorie und unter Verwendung der 
komplexen Methode werden die Durchbiegungen und Spannungen in einer diinnen 
isotropen Platte unter einer parabolischen Belastung ermittelt, die sich tiber die 
Flache einer konzentrischen Ellipse erstreckt. Der Rand der Platte ist einer elasti- 
schen Bindung unterworfen, welche Auflage und Einspannung als Sonderfalle 
enthalt. Es werden auch Grenzfalle der erhaltenen Lésung diskutiert. 
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Bemerkung zur Abelschen Integralgleichung 


Von JoacHIM FRIEDRICH, Berlin, Deutschland’) 


1. Die nachstehende Bemerkung bezieht sich auf eine Verallgemeinerung 
der Integralgleichung 


Ry , 
oe af as (1) 


yom Abelschen Typ. Gleichung (1) stellt den Zusammenhang zwischen der 


; 


seitlich beobachteten Strahldichteverteilung I(x) und der radialen Abhangig- 


_keit der Strahlungsleistung je Volumeneinheit z(7) der zylindrischen Saule 


A 


eines elektrischen Lichtbogens dar. Sie wurde zuerst von HORMANN [1]?) her- 


geleitet ; ihre Lésung i 


1) Osram-Studiengesellschaft fiir elektrische Beleuchtung. : , , 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 197. 
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ist von ihm ohne Angabe des Lésungsweges mitgeteilt und spater mehrfach 
zitiert worden [2-6]. 

In gleicher Weise ist auch der Zusammenhang zwischen der Leuchtdichte 
B(x) [cd/cm2] und der radialen Abstrahlung je Volumeneinheit J*(7) [cd/cm?*] 
eines zylindrischen Bogens durch die Gleichungen (1), (2) gegeben [7, 8]. Dieser 
Zusammenhang wird zumeist dazu benutzt, mittels eines numerischen Verfah- 
rens [1, 2, 9], vermége seitlich beobachteter Messwerte fiir [(~) bzw. B(x), aus 
(2) — oder rekursiv aus (1) — die radialen Abhangigkeiten zu berechnen. 

Im folgenden wird zunichst der Lésungsweg fiir eine verallgemeinerte 
Gleichung (1) beschrieben, anschliessend auf die Korrektur eines Ergebnisses 
fiir eine verwandte (von SupLin behandelte) Gleichung eingegangen und 
schliesslich eine Erweiterung des Wertebereichs fiir die Potenz des Kerns der 
verallgemeinerten Integralgleichung besprochen. 


2. Eine Verallgemeinerung von (1) stellt offenbar die Gleichung 


" y(t) at 
f(s) = ie (iP — she (3) 
0<a<1,p+0,a=s,a ¥ 0, f(s) stetig differenzierbar, dar. Multiplikation ~ 
mit s¢-1/(s* — z’)1-* und nachfolgende Integration ergeben 


/ PA f(s)ds  [ sP-tds y(t) dt 
waa —| wae) eae 


woraus nach Vertauschung der Integrationsgrenzen gemass einem Satz von 
DIRICHLET [10] sowie mittels 


t 


s?- ds ik bf 
/ (2? — s?)* (sP — zP)1-% > * sinan (4) 


die Gleichung 


psinan 0 [ s?-1f(s) ds 
ea : | Goa (5) 


4 


y(z) = 


folgt. Die Ausfiihrung der Differentiation nach partieller Integration liefert 
nun die Lésung 


_ psinax |, (a) Hae 
y(t) = — ae (a? — ¢P)1-% v as SF . (6) 
t 


14 

pay eae I(s)= I(x), y(t) =2ri(r), a= Ry sowie =2, «= 1/2 und 
o) = U gesetzt entsprechend Gleichung (1), dann wird Glei i 

Tat identisch mit (2). oe ‘an ae 


Ey oem 
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oy) 


3. Der verwandte Fall der Abelschen Integralgleichung 


y(t) at = 
I(s) =| (5? — Bye O<a<1;$+0, s Za) (7) 
ist von SUPLIN [11] behandelt worden. Sein Ergebnis enthalt jedoch offenbar 
einen Druckfehler; es muss richtig lauten 


p sin ax ) du 
y(@) = PER pra Le oss se Ata. (8) 


Auch in dem bekannten Werk von SCHMEIDLER [12] muss anstelle des dort 
auf Seite 214 zitierten Ergebnisses (in Beispiel 42) unsere Gleichung (8) stehen. 
Zum Beweise werde Gleichung (8) in (7) eingesetzt. Dann wird 


s s t 
_ psinex tP-1 dt fF te Vai i’ (u) du 
f(s) = ine eS oy (sP — ¢P)% (tP — qP)i-@ + | (s? — {Pye it (? — ae. ; 


Nach Vertauschen der Integrationsgrenzen im zweiten und nach Anwendung 
von Gleichung (4) in beiden Summanden ergibt sich schliesslich die Identitat 


= Ha) + | f(u) du 


Ohne Beschrankung der Allgemeinheit kann in diesem Fall a= 0 gesetzt 


_ werden. 
Mit # = 1 folgt dann aus (7) die bekannte Gleichung 


" y(t) d 
js) = [eo 9) 
0 
mit der Lésung 
t 
no = Hae (10) 4 f fee) ao 
0 ; 


gemiass Gleichung (8) [13, 10, 14}. 
Mit « = 1/2, p = 2 ergibt sich aus (7) die von v. Mises [15] betrachtete 


Gleichung 


a aa hy ) 


yty=— [1 ts ef ete reo Le foe ; (12) 


mit der Lésung 


ZAMP XI/13 
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Ein weiteres Beispiel von der Form (11), das bei der Untersuchung des 
Landestosses bei Segelflugzeugen auftritt, ist von SCHMEIDLER [12] zitiert 
worden. 


4. Es lasst sich zeigen, dass die Integralgleichung (3) mittels der Trans- 


formation 
o=a—s?, t=ar—-? (13) 


auf Gleichung (9) zuriickgefithrt werden kann. Geht vermittels (13) die Funk- 
tion f(s) tiber in g(o) und y(t) etwa in y(t), so wird aus (3) 


1 . n(t) dt 
vhs p : (@= 7)" (a? — 7)@- Dip 
und mit 
7(T) 
b (a? — r)P- Vie P(t) 
daher 


910) =f gna 
0 
Dies ist eine Gleichung von der Form (9) mit der Lésung (10), welche nun durch 
Riicksubstitution wiederum Gleichung (6) als Lésung von (3) lefert. Ein 
solcher Ansatz (13) fiir den Spezialfall = 2, « = 1/2 ist bereits von ScHMITz [2] 
verwendet worden. ; 


5. Im folgenden werde nun die Integralgleichung 


Ss 


f(s) = | (9? — ) y(d at (14) 
0 
behandelt, in der ausser dem bisher Gesagten ein Bereich (— 1 < B < n),p + 0 
und ferner eine hinreichend oft stetig differenzierbare Funktion f(s) zugelassen 
sein soll. 

Gleichung (14) verallgemeinert die bisher betrachteten Falle im Sinne einer — 
Erweiterung des Wertebereichs fiir 8; sie erfasst damit zugleich zwei weitere 
Spezialfalle. Der eine (n = 1) ist in der bereits genannten Arbeit von SuPLIN 
[11] enthalten, wahrend der andere Sonderfall (p = 1) mit einem Ergebnis von 
ROTHE [16] identisch ist, das im Hinblick auf ein Anwendungsbeispiel in der 
Wasserbautechnik entwickelt wurde. | 

Es ist bemerkenswert, dass diese Verallgemeinerung, die einen gemeinsamen 
Ausdruck fiir alle bisher bekannten Anwendungen der Abelschen Gleichung 
liefert, ebenfalls mit verhaltnismassig einfachen Mitteln zu behandeln ist. 


Vol. XI, 1960 Bemerkung zur Abelschen Integralgleichung 195 
Multiplikation der Gleichung (14) mit (z’— s?)"s¢-1 und nachfolgende 
Integration iiber s ergibt 


z 


fee ~ s0)F 50-1 f(s) a= f (20 — sb)H# sb- ft (s — 10)" y(t) dtds; (15) 


0 0 0 
hierbei sel —1 < w <n und 6+ u+1=n = 0 eine ganze Zahl. 
Wird nun das Integral auf der linken Seite von (15) etwa mit w,(z) bezeichnet 


[ (a — s*\" 51 f(s) ds = yo(2) , (16) 


6 


dann ergibt sich aus (15) durch Vertauschung der Integrationsgrenzen auf der 
rechten Seite 


pole) = EEE E ED [er — ey yt) at (17) 
0 


wegen 


fe — sh)# (sb — tP)P sel ds — 


t 


(2? os ae ae 


P 


B(B+ipti) (18) 


(wie durch die Substitution (z¢ — s¢) =u (z¢—#) leicht zu sehen ist). 
B (6 +1, + 1) ist die Eulersche Betafunktion. 


Wird ferner 


il a) 
2-1 (Og W,(2) = Wr+ai(2) (R =O, Led 1) (19) 


gesetzt, dann ergibt sich als Lésung von (14) 
1 7) 


yel= saree iy Pay on Po (20) 
oder 
= pear GENT WED 
x eo je ae (ref (0 — uP)Hab> f(a) a) ai 
Setzen wir hierin speziell p= 1, dann folgt 
Wt) = FEF TparEn mee | = tHe fe i 


ein Ergebnis von ROTHE [16]. 
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6. Es sollen auch die Differentiationen in (21) noch ausgefiihrt werden. 
Wird hierbei jetzt /(0) = 0 gesetzt, dann ist zunachst gemass (19) und (17) 


=f (to — sh)"f"(s . -"(23) 


0 


und weiterhin 


n= pe [ C4 Pe fonsorroat| cae (24) | 


Aus (24) lasst sich sofort ein weiteres Ergebnis von SuPLIN [11] ablesen, das 
hier dem speziellen Fall = 1 entspricht. Gemass (19) fiir k = 1 ist 


0 
1?) w(t) = bry Y(t) - 


Wird hierin y,(¢t) aus (24) eingesetzt, dann ergibt sich aus (20) mit u = —f 


(wegen ” = 1) die Losung der Gleichung (14) zu 


) _ Sinpx 1 a ROE. f"(u) u du R 
as at Dilger oe sa + f Rees i (25) 


fiir diese gilt also — 1 < 6 < 1 mit zweimal stetig differenzierbarer Funktion 
f(s) und keinesfalls allein 0 << 6 <1, wie von SuPLIN angegeben. Es ist zum 
Beispiel y(t) = f(¢) sehr leicht zu sehen fiir B = 0. 

Die weitere Differentiation gemass (21) ergibt nun die Lésung 


(i) be 1 1 
WT PAT BEY T+) | pee 


44 t 
«3 ton | (= wey f° 0) 
v=0 6 


(26) 


der Gleichung (14) mit den Voraussetzungen der Nr. 5. Die Koeffizienten a 
gentigen dabei den folgenden Bedingungen: 
Es ist 


ae 


a,,=\1firk=7 


(27) 


Gia, + @; pa b;, k-1 
(ki > 1,4 :0)5 


——— 
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Ferner ist 


k 
aon= [+o (u—» +1), (28) 
b,=([L+1+p(u—)]; (29) 


die Laufzahl k lauft in diesen Formeln von 1 bis n. 

Fir 2 = 1 wird gemass (28) a), = (1 + uw), gemaiss (27) wird ay, = ayy = 1. 
Damit ergibt (26) das bereits vorher abgeleitete Suplinsche Ergebnis (25). 

Ist jedoch /(0) + 0, {dann tritt auf der rechten Seite von Gleichung (26) 
noch ein zusatzlicher Summand 


Ihr NT (u+ 1) ™ (@ — 1) -«- (4 — n+ 1) eet 


auf. 

Die Untersuchungen haben gezeigt, dass auch die Lésung der hier bespro- 
chenen weiteren Verallgemeinerung der Abelschen Integralgleichung, deren 
Spezialfall eine Anwendung in der Gasentladungsphysik besitzt, mit einfachen 
Mitteln erreicht werden kann. 
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Summary 


In this paper the solution of a generalized integral equation of Abelian type is 
communicated with references to an application in gas discharge physics and, on 
the other hand, to the correction of a result in the corresponding literature. An 
extension of the validity range of this equation has been considered furthermore. 


(Eingegangen: 5. September 1959.) 
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A Linear Heat Problem with a Moving Interface 
By Rosert E. Gipson, London, England’) 


1. Introduction 


In this paper we consider a one-dimensional heat conduction problem which 
possesses some unusual features, and which may be formulated as follows: we 
seek functions u,(x, t) and u(x, 4) governed by the equations 


2 Ou, Oy . (1 
1 pve PO 0O<« < s(t); (1) 
q OU, OU 

a = —F, - st) <#eo, (2) 


and satisfying the initial and boundary conditions 


u,(0,t)= g(t), t>0 (3) 
u,(s(t), t] = u,[s(é), 4] = htt), *>0 ’ (4) 
u(x, 0) = f(x) , (5) 


where the functions g(¢), A(¢) and f(x) are given. The motion of the boundary 
x = s(t) separating the two regions is also specified, and we suppose further- 
more that s(0) = 0. The problem as stated above for the semi-infinite domain 
x > 0 consists in reality of two distinct problems, one relating to the interval 
0 < * < s() (region 1) and the other to the interval x > s(t) (region 2). Prob- 
lems of this type involving one or other region arise in a variety of physical 
situations some of which have been discussed by authors [1 — 4]?). In these — 
papers solutions in closed form are given to examples in which g(t) and h(é) are 
constants, f(x) is a linear function, and s(¢) is given by 


Sy (6) 
or 
Se Pie; (7). 
but, save in one instance [3], the techniques adopted cannot readily be extended 
to resolve problems involving boundary conditions more general than these. 


1) Imperial College, University of London. ; 
*) Numbers in brackets refer to References, page 206. } 


ee 
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We shall shew that for both regions the determination of a solution can be 
made to depend upon the possibility of solving a singular Fredholm equation 
of the first kind. In particular, when the motion of the boundary is restricted 
to the form (6) or (7) the kernel of the integral equation reduces to a type which 
enables a formal solution to be obtained in terms of a contour integral. We find 
that when the boundary moves with uniform speed the conditions for existence 
of the integral are not unduly restrictive, but when it moves according to (7) 
conditions must be imposed which are unlikely to be satisfied in many practical 
problems. For this latter case an alternative form of solution is discussed. 


2. The Temperature Distribution in Region 1 


2.1 Derivation of the Integral Equation 


In region 1 the problem reduces to the determination of the temperature 
field u(x, ¢) satisfying 


Siew ya ae O =e <2 Sif) y (8) 
a{O0, t)}=e(7), #> 90; (9) 

and 
u[s(t), 4] = h(t), t>0, (10) 


where, for convenience, the suffix has been dropped. 
It is evident from condition (9) that some simplification results by intro- 
ducing a doublet of strength g at the origin, and taking a new dependent 


variable 
t 


pau ——*_ [ 62) G (et 1) G1) ae, (11) 
2al/n F 
where 
G(x, t) = 7/2 exp (5) i 
_ Since 


u(0, t) = v(0, t) + g(t) 


_ it follows that v must satisfy 


a = ae, ORS wa SES (12) 
ase i 6 (13) 


mand 


ats), 1 = a) — SO [e(x) G[s@), #— PE ac. (14) 
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We now seek a solution of the form 
v(t, 1) = [pO (G@—-E)-GH+E DIME, (15) 
0 


where @ is an arbitrary function, and this choice ensures that (12) and (13) are 
met with. It should be noted that ast > 0 
vV>2a |x yx) , 


but this violates no imposed condition since region 1 also vanishes at this limit. 
Finally, condition (14) at the moving boundary leads to the following 
integral equation for the unknown function ¢(6): 


[oO {G ise) — £4 - GI) +80} a8 
(16) 


sae sane G(s), tJ (¢— 1dr, 


which may be written in the form 


[oe sinh [5207] exp(—zanq) B= 20. (17) 
0 


where the right side of (16) is denoted by 


27-12 exp( S ae ) x(t) . 


2.2 Solution for Region 1 when s = at 


This integral equation may be solved when the boundary moves with 
uniform speed, for then 


Sedan) eS 


and the hyperbolic function in the integrand of (17) may be absorbed in the 
unknown function. By obvious changes of variable (17) can be cast into the 
form of LapLacr’s integral equation (see [5]) with the solution 


y(&) = 2 € cosech (5) = ie (sar) exp (p &) dp, (18) | 


where the choice of the real constant c will be dictated by the form of y. This ; 
constitutes a purely formal solution to (17) under the conditions stated, and it 


eee we © 
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has been tacitly assumed that A(t) and g(t) behave in such a way that the 
contour integral (18) exists; the restrictions on g and h which ensure this are not 
severe, but we shall not pursue this point further here. 


As an example of the application of the method consider the simple case 
where 


& = (const) and A= hy (const). 


From the definition of y it follows that 


«(ap) ~ sayp Phe caryg)| oP Cseare) 08 


where the complementary error function is defined by 


The integral (18) may be evaluated by standard methods, and we quote here 
only the final result, namely 


g(é) = ——= [o coth (Sy 7) — g, cosech (74) exp (S25) ee) 


ae 


whence it follows from (20), (15) and (11) that 


co 


x 1 
w= gperie( 7) + tims ia (Far) 7) I | [Ho coth ( 45) | 


Se cosech (#°,) exp( 57) one exp(z a7) deen. | 


2.3 Solution for Region 1 when s = Bt” 
In the next case we suppose the boundary to move according to the law 
c= pre 


and in these circumstances equation (17) reduces to a form which can be solved, 
at least formally, using Mellin transforms [6]. Writing 


Ez = 
¢ 2ayt 


we have from (17) 


fi plé 2) sinh(£ C) exp(—¢*) df = (x) (22) 


0 
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where 


The Mellin transform of y is denoted by 


=f y(o) Ce ae 
0 


with a similar notation for the other variables. Thus multiplying (22) by t?~* 
and integrating over (0, co) we obtain 


#(o) = p*(0) 20-902 r (1 —o) exp(h a) [D ~B\ ye | ee eae 
uy (0) = (o) Pp 8 a2 o-1 a2 o-1 ay2 , 
where D is WEBER’S parabolic cylinder function (see for example [7]). Solving 
for y* and inverting we find 


b+1t00 
2 


Da Fie E-° de 
g(E) = 29" exp(— ue vite 


Jani | Pa=a[, ,Calay?)~ 0, Olay] 
(25) 


If, however, we consider the example of the preceding section where g and h 
are supposed constant we find, from the definition of y and (23), that 


y(t) = =e exp (qr . ax) [Ho — £% erfe( B 7) | == A Say, (26) 


the Mellin transform of which does not exist. This difficulty may be overcome 


for example by introducing an exponential factor to ensure convergence of the 
Mellin integral. We take 


y(t) = p(t) exp(—dt), O>0, 


of which the Mellin transform is 
y*(o) = AI) 6°", Re(o) > 0. 


The contour integral may now be evaluated using the calculus of residues, and - 
we find that 


A Al he lim . (—1)” gin +3)/2 5n en ‘ 
p(é) = exp(—5 8 a 4 a0 0 (n!)2[D_»_s (—B/a V2) —D_,,_, (Bla /2)] (27) 


of which only the first term contributes in the limit. The expression may be 
simplified by using the known result 


De, .4(*) = bee i (=)" exp(——) oe [ex (+) erfe(—| A 
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where 1 is zero or a positive integer, and the solution is now expressible in the 
form 


=| — 80 eric (6/2 a) 
m8) : 2 a//xerf(B/2 a) 28) 
It follows from (15) and (11) that 
rt (4/2 
4 Bo + (hy = &) EV (29) 


a result obtained previously by GRINBERG [1]. It may be noted that for this 
particular example the rather involved procedure outlined above can be’ 
avoided for the integral equation (22) may be solved by inspection. Further- 
more, if the function yp is expressible as a polynomial: 


Vy = > 4a t A, > —2, 
it is seen from (22) that 
p(E) = D7 b, E, 


where the coefficients b, are determined from 


bn 2'3+4n)/? exp (— B?/8 a”) 
a 


n I(1+4,)[D_,_;, (-B/ay2) —D_,_,, (Bla y2)] | 


3. The Temperature Distribution in Region 2 
3.1 Derivation of the Integral Equation 


In the second region x > s(¢) we seek a function u(x, ¢) governed by 


‘i 02u = Ou (30) 


ges Se TE 
which satisfies the initial condition 
u(x, Or f(x) ,. % > 0 (31) 
and the boundary conditions 
u (s(t), t) = hh), 71> 0, (32) 
St (c0, )=0, (33) 


where f(x) is assumed to be of bounded variation. 
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If a solution of (30) satisfying (33) is chosen of the form 


ene fue G(x —E,) +A GH+E, 0] aE, (34) 


then this corresponds physically to the temperature at (x, ¢) due to an initial 
temperature distribution f(x) in x > 0 and A(x) in x < 0. Since of these initial 
conditions we are concerned only to ensure that (31) is met with, we are at 
liberty to regard 4 as an arbitrary distribution which must be so adjusted that 
the remaining requirement (32) is satisfied. In this way we are led to the follow- 
ing integral equation for A(&): 


2a x h(t) = [8 G [s(t) — €, 4] + A(é) G [s() + €& ty} ae. (35) 


3.2 Solution for Region 2 when s = at 


When the boundary maintained at temperature A(t) moves into the medium 


at a uniform speed «, the above equation reduces to 


co 


2ayact Me exp fas) =| [M6 exp(zas) + AC) exe( Ga: )] e*P(aare) # 


0 


(36) 


and following again the methods of [5] it can easily be shewn that a formal — 


solution is 


As illustrative of the method we consider the case where the temperature 


at the moving boundary is invariable, namely when 
h = hy (const) . 


The evaluation of the contour integral presents no unusual features and it is 
found that 


ae) = ha [t+ exo(22)] ~ 10 ex0(42). 
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bo 
oO 
nn 


which enables the temperature distribution to be expressed in the form 
eer x h % %*—2at 
u = — erfc = > ex [-= —at ——— 
% conrad Re z &XP| —— (x — a t)| erfc 2ayi 


+ a exp( ati) [re g) sinh | * Fae] exp(par = oT) dé. 


An alternative method of approach to these problems of region 2, which 
has been exploited by BENFIELD [4] and GRINBERG [1] and which is discussed 
in [3], depends upon the introduction of a system of moving axes with origin 
embedded in the moving boundary. This is tantamount to taking a new space 
variable 


z= x—s(t), 
in terms of which the governing equation (30) becomes 


02u ds ou Ou 
0z2 dt cubs dks 


a2 


The success of the method depends upon two circumstances: the boundaries 
are now ‘fixed’ at z= 0, oo, and furthermore, since ds/dt is constant, the 
independent variable ¢ may be eliminated temporarily by means of a Laplace 
transform. It is clear, however, that no special advantage would have been 
gained by using this approach to the problems of region 1. 


3.3 Solution for Region 2 when s = B t1!” 


An argument corresponding closely to that of Section 2.3 can be appealed 
to for deriving the solution in region 2 appropriate to the case when: the 
boundary moves according to (7). By the same changes of variable as used in 
that section, the following equation for determining A may be derived 


[uc g) exp(—£ 6 = ¢*) df= pln), (39) 
0 


where 


1) = /xexp(fox) * (42) - | ecvexn( ¢-@) a, 


and using the Mellin transform as before it can be shewn that 


b+100 ; 
1 2t- 9 w*(o) exp (—B?/8 a?) &-" do 40 
Ble ree ae Tl (1.— 0) D,_, (8/42) % 


and the remarks concluding section 2.3 apply mutatis mutandis to the above 


- solution. 


ee 


(38) 
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4. Concluding Remarks 


The problems discussed in this paper may be generalized in two obvious ways. 

Firstly, the condition s(0) = 0 may not be satisfied, and in this event 
some initial temperature distribution will be prescribed in region 7 as well as 
in region 2. In region 2 a change of variable 


y = x —s(0) 


reduces the problem to that already considered above. For region 1 it would 
be natural to attempt to extend the methods used in the paper and choose as 
unknowns the initial temperature distributions outside this region, that is in 
x <0 and x > s(0). A similar problem where this procedure proved successful 
is discussed elsewhere [8]. 

Secondly, when the boundary moves in a manner other than according to 
(6) or (7) it seems unlikely that explicit solutions can be obtained by the 
methods used here since the kernel of the integral equations will no longer be 
expressible in product or difference form. An exception to this may be noted: if 


s=at+ pt 


a solution for region 2 may be found. 
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Zusammenfassung 
Die Temperaturverteilung in einem K6érper, der den Halbraum einnimmt, 
«4 = 0, wird unter folgenden Bedingungen bestimmt: 
1, Die Temperatur der Oberflache x = 0 ist eine gegebene Funktion der Zeit f, g(t). 
2. Die Temperatur einer gemadss ¥ = s(t) fortschreitenden Flache stellt eine ge- 
gebene Funktion der Zeit, h(t) dar, wobei s(0) = 0 und h(0) = ¢(0). 


33 ey = 0 ist die Temperaturverteilung des Korpers eine gegebene Funktion 
von 4%, f(x). 


£ 


Fs 
{ 
: 
j 


Die Lésung des Problems wird auf jene einer linearen Fredholm-Integralgleichung . 


zuriickgefiihrt. Es wird gezeigt, dass diese Inte 
besitzt, wenn die Funktion s(¢) proportional zu ¢ oder #12 verlauft. 
(Received: October 2, 1959.) 


gralgleichung eine explizite Lésung — 


a 
. 
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Die Harmonischen der Rotationsenergie bei der Schubkurbel 
und verwandte Fourier-Reihen 


Von WALTHER MEYER ZUR CAPELLEN, Aachen, Deutschland?) 


1. Grundlagen 


Wahrend die Koeffizienten der Fourier-Reihe fiir die Bewegung des 
Schubgliedes 3 und fiir die Winkelbewegung der Koppel 2 (Figur 1) nicht in 
Form geschlossener Ausdriicke anzugeben sind [1, 2, 3, 6, 7, 8] 2), kénnen, wie 
hier gezeigt werden soll, fiir die Rotationsenergie des Koppelgliedes, welche ja 
proportional dem Quadrat ihrer Winkelgeschwindigkeit ist, einfache Entwick- 
lungen angegeben werden. 


Figur 1 
Geschrankte Schubkurbel. 


Aus Figur 1 liest man mit, « als Kurbelwinkel, A)f = e, A4jd =a, AB =e 
und A = alc, e = e/c, wobei A + ¢ <1 sein muss, fiir den Koppelwinkel @ die 
Beziehung c sing = e + a sing oder 


sing = e+ Asina (1) 


ab. Durch Differentiation folgt hieraus fiir das Verhaltnis 


Wo dgyldt dp 
=o,  daldt ~ do 


1) Lehrstuhl fiir Getriebelehre, Rhein.-Westf. Technische Hochschule. oe . 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 217. 
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die Gleichung 


dy A cos % 
t= a esp 2) 
worin 
cos?@ = 1 — sin?’m = 1 — (e + Asina)? (3) 
ist. 


Da die Rotationsenergie des Gliedes 2 von 1/cos*@ abhangt, soll zunachst 
dieser Ausdruck analysiert werden. 


2. Vergleichsreihe 


Hierbei lasst sich die Aufgabe leicht durch Vergleich mit einer von der ~ 


zentrischen Kurbelschleife her bekannten Entwicklung [4, 5, 7, 8, 10] lésen. 
Es gilt fir 0 <7 < 15) 
7 (n + cosa) 


a = re 2 3 —_— eee = 
eet oy cose ene yy cos2. a+ 7% cos3 « ate 8 (4) 


und danach auch offensichtlich 


1 1 : 
Yo Tf yi 2ncsa  1— 7 l 21,) 
oder 
1 1 2 
1 Pim ga po icosg etl = teen eens a (— )" cosn a . (5) 


Ms DAD) cies 


3. Die Reihe fiir y = 1/cos’*p 
3.1 Partialbruchzerlegung 


Nach der von der Integralrechnung her bekannten Methode der Partial- 
bruchzerlegung lasst sich hier schreiben 


1 1 ; 
= cosig — fie Aeine)zC ee aoe (6) | 
worin 
RAL. 1 aft il 
Y1 2 l+e+Asing’ M2 = Deel eee etnies (7) 


bedeuten. Jeder der beiden Ausdriicke lasst sich jedoch auf yy nach Gleichung 
(5) zuriickfiihren. ; 


8) Es kann auch —1 <7 <0 sein, aber bei d i i 
ee cae i er praktischen Anwendung gelten die oben 


— Oe ea eee 
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3.2 Entwicklung fiir y, 


Erweitert man y, mit dem zunichst unbekannten Faktor y, und fihrt 
a = 2/2 —« ein, so wird 
fy/2 My/2 


|e i ay Ce Sp eI at > er i meee (8) 


sofern 
jy (E + 8) = 1+ of nnd pyA—2 (9) 


gesetzt wird. Addiert und subtrahiert man beide Gleichungen (9), so folgt 
durch Division 


L$ m\?. teeta 
Gases: (10) 
oder 
/j \ Beh Ree ewe 
~j Ae ES oes Bas) ee eS eee 
V1 (7 + ( )+ Y1—( — e) 
Ferner muss 
Z 7 
A= A (12) 


sein. 
Da innerhalb der gegebenen Grenzen (vgl. auch Grenzfalle Absatze 3.5c 
und 4.2c) der Parameter 7, kleiner als 1 bleibt, gilt somit 


= il v2) am 
m= lm) = laaag agp OC ntcosm a]. 13) 


Bynes 


3.3 Entwicklung fiir Vy, 


In gleicher Weise liefert die Erweiterung von y, mit #,/2 und die Einfithrung 
von a* = 7/2 + « jetzt 


ue [o/2 
Vi a Yo (Na, &*) = 1+ 7 reaee cosa* ’ ae 


und in ahnlicher Weise folgt aus der Bedingung 
fp (1L—e)=14+ 78, pbgA=2 mp (15) 
der Parameter 7 zu 


= Vi1+ @—«) —Vi—(A+ 2) a RO s— 7? (16) 
Yi+ (@—e+y1l—(+e) 


- mit 


be (17) 
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Es gilt somit 7, = 7, (—é). 
Damit hat y, die Form 


ll “a. | sl per ati oS (— 7)” cos a] (18) 


eS Yo ie =r ge ee 
n=, 


ypere 


3.4 Entwicklungen fiir y selbst 


Beachtet man nun die Umformungen 


cos2 nm 
= cos(na = 2na) = (—1)* cos2na, (19a) 
cos2 n a* 
cos(2” — 1) ‘ 
| = cos[=" : mE (2 — 1) a| = +f (-—1)*sin(2n—1)a, 
cos(2 ” — 1) a* | 
(19b) 
so fiihrt die Addition der Gleichungen (13) und (18) auf 
1 
Weta tact A, 4 ee Me, cos2 no 
4 (20) 
+ Dd) Bez-1sin(27—1)a. 
n=1,2,... 
Hierbei wird 
1 n n ih 
Ay=5 [5 ao te ed = 3 (1+); (21a) 
Ag, = (—1)" (4 13" + 2 183") » (21b) 
Beg 1)” (Wyte ree ae (21c) 
wobei die Abkiirzungen 
2 14 1 2 n il 
ek = [EAS E 2 aad 
Seay Le a V+ e2?—2’ VAP . ja -ef— mR (21d) 


benutzt wurden. 
3.5 Sonderfiille 


a) Fir A= 0 wird getrieblich die schiefe Kreuzschleife erhalten, es wird 
y = 1/(1 — e?), das heisst, es wird 
1 
Age 1— 2 


und alle anderen Koeffizienten verschwinden, wie auch fiir 4, =.= 0 aus 
den Gleichungen (21) folgt. 


EE 
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b) Fiir e = 0 liegt die zentrische Schubkurbel vor, und die Gleichungen (11) 
und (16) liefern 


me Viz A-Vl—A 1 ee 
hy ee -Y1-#) a) 
sowie 
1 
y= v= V= 1 ag (22b) 


Somit verschwinden — da ja auch y jetzt eine gerade Funktion ist — samtliche 
B, und es folgt vereinfacht aus den Gleichungen (21a) und (21b) 


1 
Ag =Vv>= yi—-# ) (23a) 


Agese (=I 29.47", (23b) 
c) Fire + A= 1 erhalt man den Grenzfall des Getriebes mit Verzweigungs- 
lage. Dann geht aber y fiir « = 2/2 nach Unendlich, die Reihe konvergiert 


nicht, wie auch aus Gleichung (21d) zu erkennen ist: Fiir e+A=1 oder 
A = 1-—e geht », nach Unendlich. Es muss also hier ¢ + 4 < 1 sein. 


3.6 Geometrische Deutung der Parameter 7; 


Die gréssten Werte des Koppelwinkels, welche fiir « = 2/2 und « = — 2/2 
erreicht werden, seien 9,,,, und 9,,;,. Fiir diese gilt 


@é—a 
Cc 


=e—/, 


SIN @max = ——_ =) ae SIN @m in = 
und setzt man dann 


aE 14 
Pax = 7 -2Y1 CunHa > 2, 


so wird 


y1 + (e+ A) = yi+ cos2 y = 2 cosy ; y1 —(e+A= y2 siny , 


Vi + (A—2) =V1—cos29=y2sind, Y1— (2-2) =//2 cose. 


Also folgt 
__ cosy — cos d sind — siny (24) 
"h =~ “cosy + cosé ’ "2 “sind + siny ’ 
worin 
I 1 ML i 
Y= a- NW ihes a 5 Pmin 
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bedeuten!). Beim zentrischen Getriebe (e = 0) wird 


Pmin = —Pmax > also P= 5 =~ 


und 
cosy —siny _ Peace 
1 = Ne 1 = “cosy + siny apes eh (25) 
Setzt man formal 
1l-+e —€E : 
= coshd, = sinho , (26a) 
A A 

so wird auch ganz einfach 

“T= 65 Nae (26b) 


A = Asinho, & = Asinho. (26c) 


V4 2 
Auch lassen sich mit der Substitution sind = A/(1 + «), siné = A/(1 — «) die 
Parameter auf die Formen 


i= tg? Ya te” es (1 + €) cos oe — =(1—e)coso (26d) 


? 
2 Vy 2 


bringen. 
Hinsichtlich der Darstellung der Parameter vgl. Figuren 2 und 35). 


4. Die Fourier-Reihe fiir w = i? 


4.1 Allgemeine Entwicklung 


Die Reihe fiir 


2 2 
ee A as eg (27) 
cos? p 


lasst sich auf die Reihe fiir y bzw. y, zuriickfiihren. Es wird 


Aa WE COs Zia 2 22 (28) 


v= 2 O cos? p = 2 (y + y cos2 a) = = (wy + w,) . 


Fiir den Anteil w, = y cos2 « hat man die Reihe fiir y gemdss Gleichung (20) 
mit cos2 « unter Beachtung der goniometrischen Formeln 


2 cos2n% cos2 «= cos2 (n+ 1) a+ cos2 (n—1)a, 
29) 
2sin(2m— 1) «sin2 «= sin(2 +1) «%+sin(2”— 3) « 


4) Diese Umformungen spielen auch bei der Bestimmung des Extremums der Winkelgeschwin- 
digkeit eine Rolle [6]. 
5) Ahnliche Tafeln kommen bei der geschrankten Kurbelschleife vor [9]. 
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Figur 2 
Funktionen 7, = f(e) und 7, = f(€) fiir verschiedene Parameter A. 


1 


°F 


Figur 3 
Funktionen 1/y, = (A) und 1/y, = f(A) fir verschiedene Parameter €. 
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zu multiplizieren und erhalt 


1 
W> = A, = (2 Ay + Ay) cos2a+ 5 (Ag + Ag) costa 
+ > (Ag + Ag) cos6 a+ ++ 
(30a) 


2 1 4 
4 = (—B, + B,) sine + > (B, + B;) sin3 « 


a S (B, + B,) sin5a+e. 


Addiert man nun die Reihe fiir w, zu der von y selbst, so erhalt man die Ent- 
wicklung 


2? cos? « 
are lt (e+ Asina)? Ay tora cos2  % 
oe (30b) 
+ 2 BF, 1 sin(2n —1)a, 
P= Ea iaeie 
worin zum Beispiel 
,2 
At = 7 [4.+2 Ag+ Ag] 
wird, das heisst 
A? ‘ 2 
At = | [— (11 nt + 2 8) + 2 (v1 1 + ¥2 3) — (r4 Hf + Ye 78), 
2 5 
A= —|y [1 mt (L — Wf)? + ve 9g (L — 18)" 
A i il i 
= —% [ ~ 98) nf + (1 — 18) nf] = —(— nt + 5 a) 
Zusammengefasst erhalt man dann 
Aj = mm (m1 + m2) 5 (31a) 
A* ait 1yr-3 (— an 2n 1b 
Bn a, amie & (31b) 
n= 1 n— 1 ' 
BE y= (=1)8=1 f= ngn—1 — ar ng"-1) : (31¢) 
Aus 7/7, = 4/2 = np/r, folgt auch 
Aga = (lps ie eee (31d). 
Vehad 1\"-1 A 2n-2 2n—2 | 
PM Gee Ye oy ale (31e) 
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insbesondere A 
ih aD Ne IE ao 
BY = (2 ze 25 ("3 — ni) -°) (31f) 


4.2 Sonderfalle 


a) A= 0 liefert trivial w = 0. 

b) «= 0, das heisst, die zentrische Schubkurbel liefert 4, = y. = y [vegl. 
Gleichung (22a)] und fiihrt auf 

AS =n, AB, = = (1-8 8, 

alle B=0. 

c) Der Grenzfall e + 4 =1 fiir das durchschlagende Getriebe fiihrt jetzt 
aber zur Konvergenz. Zunachst lasst sich w mit 
sing = (1 — Aa) + Asine = 1 —A (1 — sina) = 1 —A (1 — cosa) =1 — 2A sin? = 
auf die Form 
A cos? (%/2) A (1 + cosa) 


= T—Asin?(@/2) — (2—a) + Acosa 


(32) 


bringen. Wahrend 7 oder w = 2? gemass Gleichung (27) einen unbestimmten 
Ausdruck fiir « >a/2 liefert, zeigt jetzt w in der Form der Gleichung (32) 
keine Unbestimmtheit mehr; denn « = 2/2 oder « = 0 ergibt 


7 
w(5) — A . 
Man kénnte unter Riickgriff auf die Gleichung fiir yy hieraus die Reihe ent- 
wickeln, man kann aber auch unmittelbar e+A=1 in die vorstehenden 
Gleichungen einsetzen. 
Man erhalt zunachst nach den Gleichungen (11) und (16) 


y= es =n und y=1, 
ferner 
us iz pee ane 
also gemass Gleichung (31): 
At=1—-yl—A=1-fe, (33a) 
Ag, = (-1)-22/1— ay, (33) 
Be = (1-821 = Agate (33c) 


8) Die zweite Form ist fiir die rechnerische Auswertung ginstiger. 
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Die Koeffizienten Ax, und B¥,_, sind geometrisch abgestuft. Wiirde man wie 
in anderen Fallen [6] Kurventafeln fiir die Koeffizienten aufstellen, so lieferten 


die Gleichungen (33) die Grenzkurven. 
Im iibrigen reduzieren sich fiir A = 1 die Werte auf Aj = 1 — alle anderen 


verschwinden; denn es ist dann nach Gleichung (32) w = 1. 


5. Weitere Reihen 


Der Vollstandigkeit halber seien noch zwei weitere Reihen angegeben. 
5.1 w= —,—=Asinay 


Multipliziert man die Reihe fiir y gemass Gleichung (20) gliedweise mit 
sing unter Beachtung von 
2 sina cos2na=sin(2n+1)a—sin(2n—I1)a, 
2 sing sin(2 — 1) a= cos(2” — 2)a—cos2na, 
so erhalt man in ahnlicher Weise wie vorstehend 


Asina 


b= eine © te De fonts 2 Cae b,sin2na,  (34)% 
worin 
ni 3 
eS er, ra) er 
Bay = (1) (Ci — oe eG (35b) 
Dana (Le? (Sy 2 eae a) (35c) 
mit 
1+ 7 ae 1.— 
f= —— oder € =] es ae £ 
aL 1 aioe? 8 Tyee ee 


Fiir den Grenzfall ¢ + 4 = 1 konvergiert die Reihe nicht, wahrend fiir e = 0 


doch m = 2 und y, = 9 = €/V1 — 22 wird. Die ay, (nm = 0, 1, 2, ...) verschwin- 
den, und es bleibt 


b oa ee Sa \eaal 2 — y2@n-—1 
2n-a = (—1) ea : (36) 
SOG Acosa _ | 
‘ “cust, Sooo 
Dafir folgt entsprechend 
v= ya den Cos(2m—1)a+ b,, sin2 na, (37) 
m=, 2... WL 
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da ja 
2 cosa cos2 ma = cos(2m+ 1) a+ cos(2n—1)a, 
2 cosa sin(2 — 1) = sin2 na + sin(2 — 2a) a 
gilt. Hierbei haben die Koeffizienten die einfachen Formen 
Ganna Sage ch et), (39a) 
Ban = (—1)" (n3" — 18) . (39b) 


Fiir den anderen Grenzfall ¢ + 2 = 1 hat man keine Konvergenz, wahrend fiir 
é = 0 nur die Koeffizienten 


Agn—1 ak (—1)"-22 et (40) 
nicht verschwinden. 


6. Bestimmte Integrale 


Da fiir die Fourier-Koeffizienten der Funktion /(«) doch die Formeln 
2x 2% 27 


an Ay= | f(a) da, 2A,=| f(x) cosna da, x B, = | ja) sinn a da (41) 


ry 


gelten, kénnen die bestimmten Integrale auf den rechten Seiten fiir die oben 
behandelten Funktionen /(«) leicht angegeben werden. So ist zum Beispiel 


2 


da 
/ ek = 2 Ay (42) 


e+ Asina)? 
0 


mit A, aus Gleichung (21a). 
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Summary 


The rotational energy of a link in a mechanism is proportional to the square 
of its angular velocity. This ‘square’, in contrast to the angular velocity itself, can 
be expanded into a Fourier series, having finite values for its coefficients. The 
problem is solved by splitting it into partial fractions, leading to the harmonic 
analysis of the slotted lever mechanism. Further related series are given at the end. 


(Eingegangen: 21. August 1959.) 


Cylindrically Orthotropic Circular Plates 
with Large Deflections’) 


By JERzy Nowinsk1, Madison, Wisc., U.S.A.?) 


While the problem of large deflections of elastic isotropic plates has been 
considered in great detail, so far only three papers on large deflections of aniso- 
tropic plates have been published [1, 2, 3,]3). It seems, therefore, reasonable 
to reconsider this interesting problem using simultaneously, for the sake of 
comparison, two efficient and straightforward methods such as the pertur- 
bation method of Poincaré and a modification of the Galerkin method based 
on the energy concept. 

The first of these procedures was applied with success in the last decade 
to circular and annular isotropic plates by Chinese elasticians*) and recently 
by Nasu and CooLey [5] to elliptic isotropic plates with large deflections. 


An work has been sponsored by the United States Army under contract No. DA-11-002- 


2) Army Mathematics Research Center, University of Wisconsin. | 
3) Numbers in brackets refer to References, page 227. 


4) Cf. WeI-ZAnc CuIen [4]. For an extensive biography see the same author’s [5]. 
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The modified method of GaLerKIn has been extensively utilized by 
VoLmIr [7]. 


1. Consider a thin initially flat cylindrically orthotropic circular plate 
having the thickness / and the radius a. Let the middle plane of the plate 
coincide with the ry y-plane of a cylindrical coordinate system 7, y, z and the 
axis of elastic symmetry with the z-axis of the system. Assume that the plate 
is acted upon by a central concentrated force P and a uniformly distributed 
load gq. 

With reference to this system we have throughout the plate, in the absence 
of body force, the equations of equilibrium 


(7 Nj)’ —Ng=9, sed elihdy) 
M,—(rMySrnw'+ 2p 4 2) 
% “ . 2 2.96’ 


and the equation of compatibility of deformations 


Ie 07) 

[r (7 Ny)’ — BN, = ——$—(w')?. (3) 
Here w denotes the deflection of the plate, N, and NV, the radial and tangential 
membrane forces and M, and M,, the radial and tangential bending moments. 
R= E,/E,=y,/v,, while E,, Ey and »,, vg designate YOuNG’s moduli and 
Porsson’s ratios associated with the radial and the circumferential directions, 
respectively. As usual a prime indicates a differentiation with respect to 
a single independent variable (in the case considered, with respect to 7). 

For future references, we first note the known relations 


¥ hi ey 
M,=—D,(w" +t w/), M, = —D,(v, +), «4 


U Ny — 5p) > 


ees ( 
hE $ 
where D, and D, denote the flexural rigidities of the plate and u the radial 
component of the displacement vector. 
To the governing field equations (1), (2), (3) we must adjoin the boundary 
conditions which in a general case of elastically built-in edge and elastic 
resistance of the support appear as 


* Vv ; 
Bw=0, My +(-»,)n-0, wt (A—a)w'—0, at roan, 


tg 
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Here xf = h Eg/xy and x* = x,/D,, while x, and x, designate the elastic con- 
stants for the support with respect to a radial displacement and a flexural 
rotation. 

To economize on space, in what follows we shall confine ourselves, with no 
loss in generality, to a particular case of uniform load (i.e., P = 0) and rigid 
clamping (x, = cc) with edge fastened against radial slipping (#2 = oe). 

Assume that on account of the deformation of the plate its middle plane is 
deformed into a smooth surface (i.e., with continuously turning tangent plane). 
We assume, in addition, that the deflection as well as its first derivative are 
bounded functions throughout the region of the plate®). However, we do not 
postulate the boundedness of the second derivative and we shall discover in 
the following discussion that for k < 1, ie., for E, < E, the bending moments 
increase indefinitely with decreasing distance from the center of the plate. 
It follows that, in the case considered, the center of the plate becomes a sin- 
gular point. Such singularities brought about by the physical nature (aniso- 
tropy) of the medium are often encountered in the theory of anisotropic bodies®). 

At this point we find it expedient to transform the governing field equations 
to a dimensionless form by adopting the notation 


ne a: a _ w(0) erat ae gat 
x= ’ C= h’ = nN, = +5- N, » Ny = sy Ng» One 2hD ‘i 


(6) 


2. Let us first tackle the problem under discussion by using the method of — 
successive approximations based upon the smallness of a parameter. We 
choose for the parameter the ratio of the center deflection to the thickness, €,. 
Clearly, Cy =€,,ay in the case considered. By virtue of (4) and (6) the equations 
(1), (2), (3) reduce to 


Ny =, + 4%, , 7 (7) 
[x(x m,)’Y — Pn, + A(C’)2? = 0, (8) 
(eC)! Rs Ce eee (9) 


From now on the prime denotes, unless the contrary is stated, a differentiation 

with respect to the variable x. Moreover, 42 = 6 (h? — 72). Let us expand the 
quantities w, €, n,, mg in ascending powers of the parameter Co. We begin with — 
the second powers in ¢y in the expansions for n, and nin order to obtain the 
classical problem of linear bending of the plate for decreasing value of the 


| Clearly, w =0 at y = 0 because of the axial symmetry of the system. 
*) To illustrate this observation, we refer, e.g., to the book by LeKunitsKy [8] 


| 
: 
- 
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parameter. Hence we write 


@ = 0 Cy + Hg CF +s, M, = Be(x) C2 + gy(x) Ch +e, | 
(10) 

C= ala) Cot za(x) Bt ore, mtg tals) B+ Male) te. | 
Substitution of (10) into (8) and (9) and into the obvious boundary conditions 
w=w=u=0 at r=a, (11) 
yields a sequence of linear boundary value problems associated with terms of 
successive order in €y. We shall designate these problems as 1st, 2nd, ... bound- 
ary value problems depending whether they relate to terms of Ist, 2nd, ... 


order in ¢y. 
For the Ist boundary value problem we thus find 


2 ' 
Pp jing ieee (+4) =24, (12.1) 


2 (Obatl.. BO) <0; a(t) =< 0. (12.2) 


The solution of this system becomes 


fee 
fis 6 


(3 — k) —4 x*t+1+4 (k 4+ 1) x4] (13) 


_ and corresponds to the well-known solution of a respective linear problem for 
a clamped cylindrically orthotropic plate, cf. [8]. 
Furthermore, we reach 


hy = 4 (eel) (2+ 3) . (13.1) 

By a similar argument we can state the 2nd boundary value problem 
[x(x go)’ — RP ge + A*(zq)? = 0, (14.1) 
(1.— vg) go(1) + go(1) = 0. (14.2) 


Here the boundary condition (14.2) follows from the last condition (11) by 
virtue of the last equation (4) and the relation 


a(x) = golx) + * 8o(x) , (15) 


which is derived from (7). 
A simple calculation yields the solution of the foregoing system 


xs ees wrk 


a(x) = A, [- Bt x*-2 + WSS ELD oa array]: (16) 
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in which the notation 4, = —[4 (1 + k) a/(3 — #)]? and 


(17) 


B* = ae 


il Ties Bers aa a 2k+1—yv | 
k— 49 — k® 4 (k + 2) (k +- 1) (3k + 1) 


has been adopted. 
We now proceed to the next approximation associated with the third 
boundary value problem 


x(x 2)’ — h® 2h = 2? (x ag + Bo 24), (18.1) 
2,(0)=0, (1) = 2(1) =0. (18.2) 


The solution of this problem becomes 


sat 384 (k ae - 7) pen 2 (RBs [=e ue , f x) B 
i eatin pee is ( a 5) es Fae k : » E } (19) 
x (eres $23) re (ot 22) lh 


where 


A= (k+1)B+2(k-1)C+8D+(%+5\E+4+2(¢+1)F + 6h tee 


pes al ey RR: be uke 
8 (k + 2) (k + 5)’ 2(k+ 1) (342? +4k+4 1)? 
7 es it aoe k?— 4k — 57 
12 (k® — 49) (k® — 121)’ 64 (k + 2) (k + 4) (& + 9) (A® — 49)’ 
ee cited 3h? + 8k4+9 
8 (hk + 1) (e+ 2) (& + 3) (3% + 1) (3. k® + 20k + 25) ? 
Ge 1 


12 (b+ 1)2 (3k + 1)(2H 4 3h + 1) * 


Furthermore, we reach 


a, = 4 (9 — Rk) A* (20) 
with 

ee ee ba ; 

eee (20.1) 

x [4 B+ (k+1)C—(k-11) D+ 8E+(k+5)F4+2(k+1)G). / 


We confine ourselves to the approximation reached up to the present ' 
stage, which by virtue of (10) reduces to the cubic relation between the a 
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and the relative center deflection 


q@--, < Bs 
D,h = Z(H, Cy + a CB) - (21) 


Here «;, 7 = 1, 2, are determined by equations (13.1) and (20). Relative to the 
convergence of the procedure utilized in the foregoing discussion this problem 
has not yet been fully clarified from a systematic point of view, cf. [9]. On the 
other hand, it seems too much to require the convergence in the classical 
sense of the formal expansions supplied by the method of small parameter. 
In fact, in practical applications only the first few terms of these expansions 
are utilized and they may ultimately be even divergent without being ineffi- 
cient, cf. [10]. Thus we are led to regard the expansions concerned not as the 
power series but as series asymptotic in nature and this dispenses with an 
investigation of the convergence of the procedure. 

We remark that in the perturbation method the equations for the stresses 
may be inferred directly from the basic equations [cf. the last two equations (10) | 
without any recourse to the operation of differentiation. This fact greatly increases 
the accuracy of the results obtained by this method as compared with some 
other methods in which a differentiation of approximate expressions of the 
components of the displacement may occasionally lead to inadequate results. 

It may be now of interest to deduce a conclusion which will support our 
observation anticipated at the end of the first section relative to the singular 
properties of the center of the plate. To this purpose assume that k < 1. Then, 
in view of the appearance of terms of the (k + 1) order in x in (13), (19) and 
(24) the second derivative of the deflection with respect to x increases infinitely 
with x tending to zero. It follows that for k < 1 the bending moments at the 
center of the plate increase infinitely, cf. the first two equations (4). Moreover, 
the presence of terms of the (k — 1) order in x in equations (15) and (16) and 
of the (k + 1) order in x in equation (26) makes it evident that at the point 

x = 0 also the membrane stresses increase infinitely. 
On the other hand, for k > 1, it may be concluded by a similar argument 
that at the center of the plate both the components of the bending moment and 
the membrane stress vanish. This completes our observation concerning the 
existence of a point singularity at x = 0. 


3. We turn now to attack the problem under consideration by means of a 
_ modification of GALERKIN’s method. 

To this purpose, it is convenient to recall the equation of virtual work for 
a plate with large deflections, 


2a | (N, de, + Ny Beg + M, Ox, + My dns —9 dw) rdr=0, (22) 


0 
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in which 


1 } U 
&, = u(r) + > [w (es ie ae a n= = wry, ng — ha 


designate, respectively, the unit elongations and the curvatures associated 
with the radial and circumferential direction. Clearly, in equations (22.1) as 
well as in the succeeding equations (23) the primes denote the differentiation 
with respect to the variable 7. 
Let us now represent the membrane forces in terms of a stress function 
EF = Fi), h 
Re tae Ley Nga PT) (23) 
This representation permits us to satisfy identically the equation of internal 
equilibrium (1). We change now back to dimensionless variables and represent 
the large deflections of the plate in the form established in the linear theory 
of cylindrically orthotropic uniformly loaded and rigidly fixed circular plates 


€=[((3—k) —4x*t14 (R+1) x4]. (24) 


By virtue of the foregoing three equations the compatibility equation (3) 
leads to the differential equation for the stress function 


x2 F" 4% F" — h2 F’ = 28 (kt 1) By h® x (28@ — 28) by, (25) 


from which, by neglecting an unimportant constant, one readily obtains 
Be (ah 0! CM a8 Ag p48 Baer en Ca) ence (26) 
with 
k+1) 2 (k + 1)2 4(k +41 
pais z i$ les (k + 1) 
1 ine Ag? 8 2) pos Sia abt i 
and 
Coa s 


Cer DP 7) 


(26.2) 


x [(8%y — 56) Ap + (k +5) (vg—& —5) By + 2(k +1) 04 -2h—-1)C)], 


the last equation being deduced from the last boundary condition (11). 


We now return to the basic equation (22). In view of (4) and (22.1) this 
equation admits the representation 


1 
[|p (eer +er- Sc) -a re — 8) 5g ax | 
+ Shien aang cant 


0 


(27) 


ee 
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We identify the second integrand in brackets in the equation above as 
the equation of equilibrium (1) satisfied identically throughout the region of 
the plate in view of equations (23), We shall annihilate the remaining function 
and, simultaneously, derive an approximate value of the center deflection by 
substitution of (24) into (27). This process yields finally the important relation 


4 
Di = 2 (ah bo + oF 68), (28) 
with 

at = 4(k+1)(9—2), 

ak — 288 (k + 1) (k + 3) (03 — B) at 5 a cn aE ee | 
: re (28.1) 

+ae+4 G41) 40+ ETHEL IN GETS Bo | 

k+1 
+ e+ 4) (2h +1) (GA +5) Ca. J 


Apparently, the foregoing equation (28) corresponds to the equation (21) 
derived by the perturbation method and it may be established by inspection 
that a, = a in view of (13.1) and the first equation (28.1). 


_ 4. The developments of the preceeding sections permit us to solve an illustra- 
tive example assuming, for definiteness, that k = 2 and vy, = 1/4. By utilizing 
the relations (21) and (28) in combination with (13.1). (20) and (28.1), respec- 
tively, one reaches the following numerical relations between the dimension- 
less pressure and center deflection: 


qa 


perturbation method, Fame ve 120 0, 4-45.01 635 (29.1) 
4 
variational method, i< =7120 G + 45.11 @. (29. 2) 


The results are shown graphically in Figure 1 and exhibit good agreement. 
The difference, e.g., for g* = gy a4/D, h = 350 is less than 2%. 

The stresses in the plate are readily obtained by both methods once the 
deflection is known as well as the functions gj, 4, and F. The peak stresses 
occur on the circumference of the plate and the corresponding bending moments 
and membrane stresses in dimensionless form are shown in Figure 2 versus q*. 
- With reference to this figure, we note that thé membrane stresses indicated 
by the Galerkin method are in perfect agreement with the first approximation 
given by the perturbation method. A similar conclusion may be drawn with 
regard to the bending stresses. However, the second approximation of these 
stresses given by the perturbation method considerably changes the preliminary 
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Figure 1 


Deflection versus relative load q*. 
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Figure 2 


Moments and membrane stresses versus Ge 


ZAMP 


Vol. XI, 1960 — Cylindrically Orthotropic Circular Plates with Large Deflections 227, 


[Mla (nl nea 
26} Ale -(m) 
variational method —_// ie 
| perturbation 
22i1l method 
perturbation method, (2nd approx.) 
18+9 
| (Mma ° 
14 +7 Beh 
| , perturbation (Ist approx.) 
10+5 and variational method 
6+3 
4 t 
2 
025 050 + 700 150 
So 
Figure 3 


Radial moment and membrane stress versus relative deflection Co. 


results. For the sake of completeness, we also show the variation of the radial 
membrane stress and the radial bending moment with regard to the center 
deflection in Figure 3. 
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Zusammenfassung 


Es werden die Verschiebungs- und Spannungszustande in einer Kreisplatte mit 
grossen Durchbiegungen und zylindrischer Orthotropie mit Hilfe zweier Methoden 
untersucht, namlich der Methode des kleinen Parameters von PoINcARE und einer 


modifizierten Methode von GALERKIN. 
Ein numerisches Beispiel fiir gleichmassige Belastung der Platte, die Poissonsche 


Zahl v5 - 1/, und die Beziehung gi Uy = 4 wird gegeben. 


(Received: December 19, 1959.) 


Optimal Coordinates for a Simple Shear Flow 
over a Flat Plate’) 
By Kuo-Tat YEN, Troy, New York, U.S.A.?) 


Symbols 
X, Y = Cartesian coordinates 
% = (UgX}r), y= (Uy ¥/) 
U, V = velocity components in X and Y directions, respectively 
u == (U[Uy), v= (V/U) 
U, = aconstant velocity 
YW = stream function (V = v y) 
€,7 = parabolic coordinates 
2m = vorticity 
2W,) = free stream vorticity 
y = kinematic viscosity (u = 0 1). 
Introduction 


In Reference {1]*), KAPLUN introduced the concept of optimal coordinates 
for the boundary layer problems. The purpose was to obtain the boundary 
layer solutions in these coordinates ‘ valid in the whole flow field’. Consequently, 


a smooth transition from the boundary layer flow to the external flow can be 


achieved. 


In this paper, it is proposed to show that for a simple viscous shear of an 


incompressible fluid over a semi-infinite plate, parabolic coordinates can be 
used as the optimal coordinates. 


1 - 7 : ; 
) This research was supported in whole by the United States Air Force through the Air Force 
Office of Scientific Research of the Air Research and Development Command, under Contract 


No. AF 49(638)-23. Reproduction in whole or in part i i i 
SE aes part is permitted for any purpose of the United 


*) Rensselaer Polytechnic Institute. 
*) Numbers in brackets refer to References, page 236. 
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The approach used in this paper is somewhat different from that of [1]. 
First, it will be decided as to what are the necessary ‘conditions’ constituting 
a smooth transition from the viscous flow to the external free stream. A coor- 
dinate system will be called ‘optimal’ if a solution of the problem can be found 
to satisfy the necessary conditions for the smooth transition as boundary con- 
ditions. It is implied that this solution of the Navier-Stokes equations, which 
are not invariant with respect to a (conformal) transformation of coordinates, 
will not necessarily be exact, and different solutions may be obtained when the 
method of the approximation used in obtaining the solution is applied to the 
Navier-Stokes equations in different coordinate systems. 

For example, for the boundary layer flow of a uniform stream over a flat 
plate the solution in parabolic coordinates obtained by CARRIER and LIN [2] 
can be shown to satisfy the condition that the transverse velocity component 
vanishes at infinity and may be called ‘optimal’. For a more detailed discussion 
of this problem reference may be made to [3]. 

This consideration can be extended to a shear viscous flow over a flat plate. 
When the shear flow is of the boundary layer type, the displacement effect of 
the boundary layer may be important because it appears to influence the 
boundary conditions at the edge of the boundary layer. The question arises as 
to what are the proper boundary conditions to use. Thus,.a study of the 
optimal coordinates appears to be desirable. 


Coordinate Transformation and Basic Equations 


Let X and Y be the Cartesian coordinates with the leading edge of the flat 
plate as the origin and the plane of the flat plate as the X-axis. The plate is 
supposed to be situated in a simple shear flow with constant vorticity. For the 
present problem, it is convenient to reduce the physical quantities to non- 
dimensional form in the following manner. Let 


Ue Xe foe Y au) 
Mig a5 
and 
ve, 9) =— Pw, 9), (2) 


where W is a stream function defined such that 


ie Wi Wine ty (3) 


_ The non-dimensional velocity components are 
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The Navier-Stokes equations can be expressed by 
Wy Ke WY, — Yx V? Wy a Vp ° (5) 


Away from the plate, the free stream conditions are 


uw=1+2a@), (6) 
o=0, (7) 
De ow (8) 


where w, is a non-dimensional vorticity number. The conditions for smooth 
transition of the flow may be taken as (6), (7) and (8). However, it will be seen 
later that the condition (7) has to be relaxed. The boundary conditions at the 
plate are 

“=V=0, V=0- (9) 


Let the coordinate transformation be given by 
a= x(6,n), y= (6,7) - (10) 


The velocity components w and v in the Cartesian coordinates can be expressed 
as [4]: 


: 

Umeha ta Vi etae ee 
1 

hy g (= Vy Pet+ VEY) > (12) 


where J = x,y, — x, ye. The differential equation for the stream function w 
in the new variables € and can be obtained from equation (5). 
Parabolic coordinates are defined by 


e= Ei — 72, ya Qin. (13) 
Equation (5) becomes 
Yy Me — Ye WO, = Wee ap Diy » (14) 
where w is 
ae. 1 
Wee ey (Pee + Poy) (15) 


In the parabolic coordinates, the flat plate is at 7 = 0. The free stream condi- 


tions and boundary conditions can be expressed in terms of the parabolic 
coordinates by the above relations. 


— 
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Construction of Solutions in Parabolic Coordinates 


Consider first the behavior of the solution at large 7. In view of (6), the 
solution at large 7 may be written in the following form 


YP = 2 wo & fo(y) + Aly) + fol) +-°-, (16) 
and, as 7 > oo, the functions /, and /, become 
fly) 24? — Bsn — Br, hin) +2y-8,, 


where f,, P. and Pf, represent the displacement effects. Since the above form of 
the solution does not predicate to be the general one valid over the whole 
field, it is not expected that /,, B, and fh, are constant over the whole field. 
In fact, their values may change from region to region for which different 
methods are used to find the solutions f, and fy. 

From equations (11) and (12), it is found that, for large 7 


1 i 
M14 40989 —-eroe fy Bint 2Prmoen t+ Bemoe + 2m, (17) 


s 


D> aera [FPL E + 2 Ba tod? + By evo 8 7} (18) 
Thus, for << y, the velocity component u approaches 1 + 4 wy € 7 — 2 Bz wy §& 
as ‘7 approaches infinity. If the parabolic coordinates are the optimal ones, the 
solution obtained in these coordinates should yield 6,=0, because the displace- 
ment effects should vanish far upstream of the plate. This can be shown as 
follows. For § << y, the stream function y can be expressed in the following 
form 


y = Eh(y) + 208? 4 F(n) + & fa(n) ++, (19) 
where, as 7 approaches infinity, 


h(y) >2y—- fi, (20) 
F,(y) > 2 — Bs- (21) 


Substitution of y into equation (14) yields, when all higher order terms in & 
are omitted, 
fh (f—-1&)— 27k’ =0. (22) 


When 7 approaches infinity, F’ > 0. It follows that 6; = 0. Hence, far up- 
stream, the smooth transition of the flow is achieved. 
: For large é, the smooth transition is possible only when Bar= By = OAs 
will be seen later, this is found to be true. Since v approaches zero and w 
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approaches w, as 7 approaches infinity, it might be concluded that all the 
conditions (6), (7) and (8) are satisfied. However, for § ~ 7, the expression for 
v contains the term fy w, £2 4/ (£2 + 7?) which is singular as 7 > oo and hence 
is incompatible with smooth transition. This point will be analyzed later in this 
paper. 

From the above consideration, similar solutions can be constructed in the 
following manner. Let the similarity variable ¢ be defined by 


pome Dh 
6¢= = 


Ek 
and the solution be of the form 


p= 2 eri) oo aba oe (23) 


where f, and /, are required to behave as fy > 4, f; > 2, when ¢ approaches 
infinity. Then, for large 7, the above solution for y agrees with that given in 
equation (16). The index k is to be determined from the condition that the 
viscous terms and inertia terms in the Navier-Stokes equations have the same 
order of magnitude. 

Now two forms of asymptotic solutions for large € can be obtained according 
to the stipulation whether the first or the second term in equation (23) is 
considered as the dominating term. For small values of w, € the second term 
in equation (23) may be taken as the first approximation of the series expansion 
solution. Thus 


y = § fly) + 2 9 & f(y) + & fs(77) + +++, (24) 


since it is found that & = 0. Substitution of (24) into equation (14) and equating 
equal powers of € yield the following equations 


PEK + ht 0, 25) 
R+hh -hht+2ih=0, (26) 

etc. The boundary conditions are 
Lie Ae Booey (27) 

h=h=0, =0. (28) 

+4, 900, : (29) 

h=h=0, n=0, (30) 

: 


etc. Equation (25) under the boundary conditions (27) and (28) yields the 


known boundary layer solution in parabolic coordinates for a uniform flow. 
over a flat plate [2]. : 


tt tee ne 
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On the other hand, for large values of &, an asymptotic series solution may 
be attempted. It is found that, for this case, k = — 1/3, and the solution is of 
the form 


Y = 2 adq $49 ga(C) + E79 gy(C) + go(l) +-->. (31) 
The function g, satisfies 
8 @ mn m 
3 8282 +e" =0, (32) 


b> 4, Cie, (33) 
f= =0, C=0. (34) 


The solution for g, is 2 ¢?, which represents a simple shear flow. Thus, since 
Bs = Bs = 0, smooth transition of the flow is also achieved at large € The 
function g, satisfies 


3 


Hog 8 +2828 + Bb +5 es Bi + oo = 9, (35) 

and 
eet fons (36) 
=e =90, C=0. (37) 


Detailed analysis of the above solutions will not be given here, because the 
solutions of equations (26) and (35) under respective boundary conditions can 
be expressed in terms of the known ones in the Cartesian coordinates [5] [6]. 
By the transformation 7 = 7*/2, equation (26) becomes 


2h +hh—-hk+2hh=0, (38) 


with the boundary conditions fj >1, 4* > co; f,=f,=0, n* =0. The 
solution for /, satisfying the above conditions has been obtained by GLAUERT(5]. 
The value of f, is 1-73, while that of 6; as estimated from GLAUERT’s results 
is about 3.4 (GLAUERT argued that f, = 2 6, = 3-46 in our notation). 

In an analogous manner, let 


a= (za) Sle) (3) 
Coy laeaig) on 6? a 


_ Equation (35) becomes, in view of the fact that g. = 2 ee 
Gu de #2 Gt fe ce G{ = Ge 


234 Kuo-Tat YEN ZAMP 


or 
Gil + OG! —C* G+ G,=0. io 


The boundary conditions for G, are Gj >1, ¢* > 00, and G; =G)=O0rat 
¢* — (0. The solution for G, has been analyzed in [6]. It may be added, however, 
that the asymptotic series solution does not appear to be valid near the leading 
edge of the plate, while the series expansion solution becomes inadequate far 
downstream. Thus, to obtain a complete solution of the problem, a matching 
or correlating process is used in [6]. 


A Modified Coordinate System 


As shown above, the parabolic coordinates can be used to obtain solutions of 
the shear flow such that a smooth transition of the shear flow to the free stream 
flow is achieved. However, for the region near the leading edge of the plate, 
the displacement effects introduce the contribution 3 w, & y/(&% + 7?) in the 
velocity component v. This term becomes singular as 7 oo, when § ~7, 
since [3 is expected to be close to the value of 3.4 as obtained from GLAUERT’S 
work. Therefore, it appears that a modification of the coordinate system should 
be made for § ~ 7, but not at large € or 7, to remove this singular behavior of 


v. For simplicity, consideration will be limited to € > 7. A coordinate system ! 


which meets these conditions and, in addition, yields y = 0 when 7 = 0 is 


x= EF — 7, (42) 


y=2Ent+ > Ban: (43) 


For large x, ¢ is essentially «'/" and 7 is essentially yy/(x1/? + B,/2). It is noted 
that this coordinate system is no longer orthogonal as a simple orthogonal 
coordinate system satisfying the necessary conditions has not been found. The 
functions 7; and f, in the expression (24) for the stream function w are found 


to satisfy the same respective differential equations and boundary conditions 


(25)—(30). Hence in equation (43), By is 3.4. As ” approaches infinity, the 
velocity components u and v are 


1 
A(t Pope 
x {2 Bi + 8 Br E+ Prt 4 Bog (E+ En +272) +--4 


1 
"4 (P+ 0) + Bye 


u>1-40,F 


(44) 


, 


v 


j1 (45) 
x [2 + 2 BLE + 8 Be cod? — By + 28s wo F (Bs + Bo) + ---4. 


a a LA 
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The skin friction distribution on the plate can be obtained from equations 
(11) and (24), and is given by 


2 1/2 
Pon vee 


: 1 
4 (1 + 2 Bs) 


As seen from equation (43), the new y coordinate is stretched compared to that 
given by the parabolic coordinates. Thus, the skin friction on the plate is 
reduced accordingly. 

It is evident that the ‘free stream conditions’ for the region near the 
leading edge of the plate has direct bearing on the skin friction distribution on 
the plate. It is also noted that near the leading edge of the plate, v approaches 
approximately the value of /, @» + (83 @ /4). Therefore, the condition (7) for 
smooth transition is actually not satisfied. It may be possible to modify 
further the coordinate system in order to satisfy this condition, but the skin 
friction distribution on the plate will certainly be changed. Thus, it appears 
reasonable that different conditions for smooth transition for the region near 
the leading edge of the plate may be associated with physically different 
problems. Therefore, for the region near the leading edge of the plate the 
condition for smooth transition for v may be modified so that v pp aa 
a constant value as 7 approaches infinity. 


x {f(0) + 29 xt 2 (0) +++}. (46) 


Some General Remarks 


The problem of the boundary layer on a flat plate in a stream of constant 
vorticity has been analyzed by GLAveERT [5] and Li [7] by employing the 
Cartesian coordinates. In particular, GLAVERT proposed to use the boundary 
condition 

“=stl+20,9 7 =O, (47) 


‘instead of 


uw >1l+2ayv, yoo. (48) 


Of course, when the solution of the problem is expressed in the variable 
n = Y (U)/v X)1?, there is no difference between the two conditions (47) ad 
(48). The solution would be defined only for positive x. 

By using the parabolic coordinates, it has been shown that the solution 
can be extended upstream of the leading edge of the plate such that the free 
stream condition u=1-+2,¥y is satisfied far upstream. By the analysis 
given earlier for the differential equations (26) and (35), it is seen that the 
solution in. Cartesian coordinates obtained by GLAUERT and that in parabolic 
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coordinates yield the same skin friction distribution on the plate: 


Tob (eee ia {332 + 1.59 Wp ( 


mx" 5. 


y Wi: 


According to equations (17) and (18), the condition (47), although mathemati- 
cally convenient, may not be quite realistic physically. In fact, (48) appears to 
be more satisfactory, provided displacement effects are to be admitted and 
one does not interpret equation (48) to mean that uw shall be constant at 
fixed y. The fact that the skin friction distribution on the plate depends 
intimately on the free stream conditions for the region near the leading edge 
of the plate demonstrates the difficulty of imposing proper boundary or initial 
conditions at x = 0. 
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Zusammenfassung 


In eine ebene inkompressible Strémung mit konstantem Geschwindigkeitsgra- 


dienten quer zur Strémung sei eine ebene Platte parallel gelegt. Als giinstiges : 
Koordinatensystem werden parabolische Koordinaten betrachtet und damit 4hn- — 


liche Lésungen gefunden. Der Ubergang von der Grenzschicht zur dausseren Stré- _ 


mung vollzieht sich glatt, und es kénnen Wandschubspannungen und Verdran- © 


gungseffekte bestimmt werden. Vergleiche mit anderen Lésungen, die mit kartesi- 
schen Koordinaten gefunden wurden, werden angestellt. 


(Received: October 26, 1959.) 
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Zur Eiskeimbildungsfahigkeit von Kupferoxyden und Kupfersulfiden 


Von Urricn Katz, Ziirich?) 


1. Einleituns 


Im Zuge der systematischen Ermittlung der Aktivitat der verschiedensten 
Substanzen als Eiskeimbildungskerne gilt den Kupferoxyden und Kupfersulfiden 
spezielles Interesse, da sich schon bei friiheren Messungen der Einsatztemperatur 
der Wirksamkeit von Cu,O und CuS (PRUPPACHER und SANGER [8]2), Fukuta [4]) 
die verhaltnismassig hohen Werte von —4 bis —5°C ergaben. Ferner bestanden 
Zweifel an der oberflachlichen Reinheit einiger dabei verwendeter, auf nassem Wege 
erzeugter Substanzen; bei Cu,O ware vor allem damit zu rechnen, dass die Ober- 
flache schon zu CuO oxydiert ist und eine Aktivitat von Cu,O als Eiskeimbildungs- 
kern nicht oder nur teilweise erfasst wird*); ebenso kénnte an Stelle von CuS bei 
Herstellung auf chemischem Wege ein Gemisch von beispielsweise CuS, Cu,S und 
S entstehen. Deshalb schien es angezeigt, die erwahnten Stoffe direkt: in Form 
reinster makroskopischer Kristallite herzustellen. Juha 

Wahrend bis anhin im allgemeinen als Mass der Eiskeimbildungsaktivitat einer 
Substanz nur deren Einsatztemperatur gemessen wurde, ist es vor allem auch im 
Hinblick auf eine Anwendung zum Zwecke der Wetterbeeinflussung wesentlich, den 
Anteil eines kiinstlichen Aerosols, welcher zur Eiskristallbildung Anlass gibt; als 
Funktion der Temperatur zu kennen; auch iiber den Mechanismus der Eiskeim- 
bildung kénnen auf diese Weise eher Aufschliisse gewonnen werden als mit blossen 
Schwellentemperatur-Angaben. Allerdings wird fiir AgI schon seit langerer Zeit 
meist die Anzahl aktiver Teilchen pro Gramm Substanz in Funktion der Temperatur 
angegeben [9, 11, 12]; neuestens findet man auch bei SERPoLAy [10] fiir: einige 
Oxyde das Verhiltnis wirksamer Teilchen zu deren Gesamtzahl. 


2. Eigenschaften und Herstellung der Substanzen‘) 
2.1 Cu,O und CuO 


Das roteKupferoxydul, das dem Mineral Cuprit entspricht, bildet einen eigenen 
_kubischen Gittertyp, wahrend CuO, das schwarze Kupferoxyd (Tenorit) in einem 
-triklin deformierten Steinsalzgitter kristallisiert. Die Cu,O-Oberflache adsorbiert 
_zuerst O, aus der Luft, spater tritt eine sehr langsame oberflachliche Umwandlung 
in das CuO-Gitter ein [5]. Aber auch sehr lange an Luft gelagertes Cu,O liess nach 
den Beobachtungen von BrwteEscu [2] mittels Elektronenbeugung (cen feststell- 
‘bare Schichtdicke: etwa 70 A) keine Oxydation erkennen. 


1) Laboratorium fiir Atmospharenphysik an der ETH. 

2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 244. 

) 3) Dieser Gedanke wurde auch in einer privaten Mitteilung von J. H. Scuurman (Columbia 
- University) gegeniiber dem Leiter des Laboratoriums fiir Atmospharenphysik vertreten. 

| 4) Unter Mitwirkung von M. BAup (im Rahmen seiner Diplomarbeit). 


238 Kurze Mitteilungen — Brief Reports -- Communications bréves ZAMP 


Da CuO oberhalb 1110°C an der Luft zu Cu,O dissoziiert, ware es méglich, 
massives Cu,O durch Schmelzen (1230°C) des im Handel erhaltlichen Pulvers zu 
bilden, wenn nicht alle Tiegelmaterialien dabei zerstért wiirden. Das Cu-O-Zu- 
standsdiagramm [5] zeigt, dass aus metallischem Kupfer direkt durch Erhitzen an 
der Luft im Temperaturintervall von 1026 bis 1065°C das Oxydul entsteht ; als 
giinstiger Wert erwies sich 1030°C, wobei 0,8 mm Kupferblech in 5 Stunden ganz- 
lich durchoxydiert war. Beim Erkalten bildete sich an der Oberflache CuO, weshalb 
die Abkiihlung im Vakuum zu erfolgen hat. Die Priifung dieser Cu,O-Proben 
mittels Schliff, Rontgen- und Elektronendiffraktographie zeigte absolut homogene, 
bis 1 mm grosse Kristallite ohne Spuren von CuO. 

Die Eigenschaften des im Handel erhaltlichen Kupferoxyds (pro anal., Merck 
& Co.) gaben keine Veranlassung, dieses selbst herzustellen, da eine weitere Ver- 
anderung an der Luft nicht eintritt. 


2.2 Cu,S, CuyS, und CuS 


Diese unter Normalbedingungen stabilen Verbindungen entsprechen den fol- 
genden natiirlichen Mineralien: 


B-Cu,S Chalkosin, orthorhombisch 
Cu,S;  Digenit, kubisch 
CuS Covellin, hexagonal 


Weitere, nur auf nassem Wege erhaltliche Polysulfide sollen hier nicht zur Diskus- — 
sion stehen. 

Bei allen drei Substanzen weist die Léslichkeit den sehr geringen Wert von 
ungefahr 10-14 Mol/l auf. Wahrend CuS an der Luft ausserst langsam zu Sulfat 
oxydiert, sind Cu,S und Cu,S, in der Atmosphare sehr stabil [5]. 

Zur Herstellung wurden als Ausgangsmaterialien Elektrolytkupfer-Pulver und 
Schwefelbliiten verwendet; das Kupfer reduzierte man im H,-Strom bei 400°C und 
der Schwefel wurde bei 100°C im Vakuum getrocknet und von SO, und H,S 
befreit. Fir Cu,S und Cu,S, erwies es sich als notwendig, die beiden Edukte in 
st6échiometrischem Verhdaltnis innig zu mischen und bei einem Druck von 1000 at 
zu Pillen zu pressen, um diese sodann in einem Quarzrohr unter Vorvakuum ein- 
zuschmelzen. Anschliessend erfolgte die Behandlung im Ofen: in 60 Minuten auf- 
heizen auf 1200°C, welche Temperatur wahrend 20 Minuten beibehalten wird; 
hierauf abkiihlen in einer Stunde. Die Réntgenanalysen ergaben gute Uberein-— 
stimmung mit den zuverlassigsten Literaturwerten [13]. 

Zur Darstellung von CuS geniigte es, Kupfer und Schwefel im Uberschuss je in 
ein Quarzschiffchen zu bringen und diese hintereinander in ein Quarzrohr zu legen, 
welches nach Evakuieren auf 10? Torr zugeschmolzen und wahrend 12 Stunden bei 
450°C belassen wurde. Wie auf Grund des Zustandsdiagramms erwartet, war das — 
Cu-Pulver vollstandig durchschwefelt worden und das tiefblaue Pulver ergab die 
rontgenographischen CuS-Werte [1]. 


3. Messmethode und Apparatur 
. es - . : 
Die Giite eines gewissen Aerosols als Eiskeimbildungskern bei einer Temperatur 
T sei q(T) und bedeute den Quotienten aus der Anzahl Kernteilchen, die zur Bil-— 
dung eines Eiskristallchens Anlass geben, und der Gesamtzahl der Partikeln, welche 
in einem unterkiihlten Wassernebel der Temperatur T eingeimpft werden. 


—_ so pe “ots 
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In der vorliegenden Arbeit ermittelt man q(T) fiir die oben beschriebenen Sub- 
stanzen direkt in der Definition entsprechenden Weise, indem in einer Kiihlkammer 
eine Wolke erzeugt, eine bekannte Anzahl Teilchen in diese hineingespritzt und 
hierauf die Zahl der sich bildenden Eisteilchen bestimmt wird. 

Der zentrale Teil der hiezu beniitzten Apparatur (Figur 1) ist die Wolken- 
kiihlkammer, die urspriinglich lediglich zur Messung von Schwellentemperaturen 


I 


Figur 1 


Schema der Wolkenkammer. 1 Isolation; 2 Verdampferrohr; 3 zirkulierende Kihlluft; 4 Ventilator 

fiir 3; 5 isolierende Verglasung; 6 Beobachtungsrichtung; 7 Thermoelemente; 8 Sedimentations- 

zelle; 9 Impfrohr; 10 Luftauslass; 17 Beleuchtung; 12 Plexiglasdeckel; 73 Nebeleinlass; 14 Aerosol- 

Einspritzung; 75 zum Kompensationsschreiber; 76 Pressluft; 17 Eiswasser; 18 Verneblerdiise; 
19 destilliertes Wasser. 


gebaut worden war; der eigentliche Versuchsraum (Volumen 1001) wird auf drei 
Seiten und am Boden von gekiihlter Luft umspiilt, was einen kleinen, aber nicht 
verschwindenden vertikalen Temperaturgradienten ergibt (0,02 bis 0,1°C/cm je 
nach Temperatur). Zur visuellen Beobachtung der Vorgaénge im Innern ist die 
vierte Seite isolierend verglast. Ein dreifacher abnehmbarer Plexiglasdeckel ermég- 
licht die Beleuchtung senkrecht zur Beobachtungsrichtung durch ein 8 cm breites 
Strahlenbiindel einer horizontal verschiebbaren Hg-Hochdrucklampe. 

Die Erzeugung der unterkiihlten Wolke erfolgt durch mechanisches Zerstauben 
destillierten Wassers an einer mit gekiihlter und elektrostatisch sehr wirksam gefil- 
terter (vgl. List [7]) Pressluft beschickten Diise; in einem vertikalen Zylinder wird 
der Nebel auf 0°C gekiihlt und von zu grossen Tropfen befreit, worauf er in die 
- Kammer gelangt, aus der er die vom vorhergehenden Versuch verunreinigte Luft 
verdrangt. Nach der Einfiihrung des Nebels verstreichen etwa 90 Sekunden, bis sich 
die Turbulenz gelegt hat und die normale Temperaturverteilung wiederhergestellt 
ist; da die Trépfchen teils absinken, teils verdampfen, bleibt die Wassersattigung 
nur noch wahrend rund vier Minuten erhalten, in welcher Zeit die Messung durch- 
gefiihrt werden muss. 
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Ein Kompensationsschreiber registriert die von zehn in der Kammer verteilten 
Thermoelementen gemessenen Temperaturen. 

Die dosierte Eingabe von Kernpartikeln in die Wolke geschieht am besten in 
Luft suspendiert. Bei vielen Substanzen in Pulverform neigen die Partikelchen zu 
Konglomeratbildung, was die Dispersion in Luft erschwert. Mit der hier verwen- 
deten Vorrichtung (Figur 2) konnten dennoch gute Resultate erzielt werden. Zwei 
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Figur 2 
Schema der Impfvorrichtung. J Zufuhr reiner Luft; 2 Kihlfalle zur Trocknung; 3 Durchfluss 
messer; 4 Hahnen; 5 Dosiereinrichtung; 6 Zerstaubergehduse; 7 rotierende Biirsten; § Weg der 
Staubteilchen; 9 Luftstrom; 10 Glaszylinder; 11 Filter; 72 Zufiihrungsrohr in die Kammer; 13 Plexi- 


glasdeckel; 14 Sedimentationszelle; 15 Offnung zum Impfen; 16 Luftabfuhr; 17 Zapfen, der in die © 


Kammer hineinragt. 


zylindrische Nylonbiirsten rotieren um ihre horizontale Achse in einem entspre- — 


chenden Gehause, das nur gegen oben gedffnet ist und auch das zu dispergierende 


Pulver enthalt; wenn die Borsten in horizontaler Lage sind, treffen sie auf ein — 


Hindernis, worauf ihre Enden stark beschleunigt und die daran haftenden Partikeln 
in den vertikalen Luftstrom geschleudert werden, der seiner Geschwindigkeit 
gemass nur Teilchen unterhalb eines bestimmten Grenzdurchmessers mitnimmt. 
Die hier zur Verwendung gelangende Luft wurde bereits elektrostatisch gefiltert 
und in Kihlfallen getrocknet. Aus dem Aerosolvorrat im Zylinder iiber dem Pulver- 
zerstauber wird ein vertikales Rohr (Durchmesser 16 mm) gefiillt, das von oben ins 
Innere der Kammer fiihrt. Mit einer Dosiereinrichtung kann nun wahlweise aus 


diesem Rohr ein bestimmtes Volumen in die Wolke gespritzt oder in eine Sedimen-_ 


tationszelle abgefiillt werden. Letztere besteht aus einem Metallzylinder (@ 14 mm) 


' 


; 


; 
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mit verschiebbarem Kolben; ein Objekttrager (Glasplittchen oder Formvar-Haut) 
am Boden der Zelle fangt die darin ausfallenden Partikeln auf. Um die Konzen- 
tration des in die Wolke geblasenen Aerosols zu kennen, muss ein reprasentatives 
Luftvolumen aus dem Rohr in die Sedimentationszelle gefiillt werden: in praxi 
wird vor und nach dem Impfen der Wolke je eine Probe entnommen und dann 
geeignet interpoliert. Die Einfiihrung des Aerosols in den Nebel erfolgt im oberen 
Drittel der Kammer durch horizontales Einblasen von 104 bis 107 Partikeln in Ons 
bis 20 cm* Luft; da diese schon im vertikalen Rohr ungefahr Kammertemperatur 
angenommen hat, wandert das eingespritzte Volumen in Form eines Wirbelringes 
horizontal gegen die Mitte der Kammer, wobei dieser und die daraus herausfallenden 
Eiskristallchen im parallelen Licht gut beobachtet und gezahlt (bzw. geschatzt) 
werden kénnen. 


4. Messungen und Diskussion der Resultate 


Die Substanzen wurden in der Achatschale von Hand zu Pulver zerrieben; da 
die Messung an einem Stoff etwa 3 Tage beanspruchte, erfolgte das Zerkleinern 
10 Tage vor Messbeginn, damit sich eventuell an der Oberflache abspielende Ver- 
anderungen nicht im Laufe der Messung st6rend auswirkten. 

Bei der mikroskopischen Auszaéhlung der Sedimentationsproben wurden die 
Gréssenverteilungen ermittelt, die in Figur 3 dargestellt sind; kleinere Schwankun- 
gen waren dabei nicht zu vermeiden, doch ergab die stichprobenweise elektronen- 
mikroskopische Untersuchung keine Teilchen mit Durchmessern unterhalb etwa 
0,5 uw. Diese Begrenzung erschien uns sehr wiinschenswert, da einerseits die Sedi- 
mentationsproben rasch hergestellt sind, und andererseits, weil auf Grund theore- 
tischer Uberlegungen (z. B. FLETCHER [3]) unterhalb des erwahnten Grenzdurch- 
messers eine immer starker werdende Abhangigkeit der Aktivitat von der Partikel- 
grdsse zu erwarten ist. 

Beim Einblasen der Impfstoffe in die Trépfchenwolke zeigte es sich, dass nach 
2 bis 5 Sekunden die Eiskristallbildung am intensivsten und nach etwa 10 Sekunden 
bereits auf ein vernachlassigbares Mass abgesunken war. Ferner konnte keine 
wesentliche Abhangigkeit der Aktivitat von Radius und Konzentration der Nebel- 
trépfchen beobachtet werden, solange deren Vorhandensein 100% relative Feuchtig- 
keit gewahrleistete; nach der Auflésung der Wolke jedoch traten Eiskristalle nur 
noch in sehr reduziertem und verzégertem Masse auf. Die vorliegenden Resultate 
beziehen sich aber ausschliesslich auf den Fall der Wasserdampfsattigung. 

Figur 4 und 5 zeigen g(T) logarithmisch aufgetragen; die starke Streuung der 
einzelnen Messpunkte riihrt vor allem von der visuellen Abschatzung der Anzahl 
sich bildender Eiskristalle her. Infolge nicht zu vermeidender Turbulenz in der 
Kammer muss mit einer Unsicherheit der Temperaturmessung von max. + Osa Gi 
gerechnet werden. Aus apparativen Griinden konnten g-Werte unterhalb 10° nicht 
mehr gemessen werden, doch lassen die Kurven das Erreichen einer oberen Grenz- 
temperatur (Schwellentemperatur) erkennen und abschatzen. Fiir Cu,O und CuO 
(Figur 4) ergaben sich praktisch die gleichen Schwellentemperaturen, wie sie schon 
PRUPPACHER und SANGER [8] fanden, wahrend Fuxvuta [4] fiir beide Stoffe etwa 
—4°C angibt. Ausser bei —35°C, wo beide Substanzen die gleiche Aktivitat zeigen, 
lasst der deutliche Unterschied den Schluss zu, dass bei Cu,O die Partikeloberflache 
noch nicht wesentlich oxydiert war. Bei den Kupfersulfiden (Figur 5) fallt die 
nahezu gleiche Schwellentemperatur auf, wahrend sich die g-Werte zum Beispiel 
bei —10°C um Gréssenordnungen unterscheiden. Die Grenztemperatur von CuS 
differiert mit —5°C nur leicht von den Werten anderer Autoren (—4°C [8], —7,3 [4]). 


ZAMP XI/16 
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Figur 3 


Relative Haufigkeit der Teilchen in Funktion ihres Durchmessers D, (Die Ordinate ist in willktirliche 
Einheiten geteilt.) 


Infolge des einsetzenden spontanen Gefrierens von Nebeltrépfchen unterhalb 
von etwa —33°C konnten die Kurven nicht bis —40°C ermittelt werden, womit 


auch tiber den Anstieg der Wirksamkeit unterhalb —30°C keine Aussagen zu 


machen sind, insbesondere nicht iiber die interessante Frage, ob bei CuS und Cu,O 
ein Grenzwert der Aktivitat bei etwa —30°C erreicht wird. 

Fiir AgI-Rauch haben verschiedene Autoren, wie SANo [9], SmirH und HEFFER- 
NAN [11] und SouLaGE [12], eine zur unseren analoge Kurve angegeben, und FLE- 
TCHER [3] interpretiert den Anstieg der Aktivitat mit sinkender Temperatur als 
Folge der Gréssenverteilung der Teilchen; wenn auch eine starke Temperaturabhan- 
gigkeit der Aktivitaét bei Teilchengréssen unterhalb 0,1 « durchaus verstandlich ist 
und experimentell von HosLEr [6] und Sano [9] bestatigt wurde, so kénnen doch 
FLETCHERS Berechnungen in unserem Falle keine hinreichende Erklarung fiir den 
Anstieg iiber ein Intervall von mehr als 15°C liefern, sondern es ware danach nur 


t 
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CuO in Funktion der Temperatur T. und CuS in Funktion der Temperatur T. 


ein solches von etwa 4°C zu erwarten; dies um so mehr, als hier kristalline, eckige 
Partikeln mit ebenen Flachen zur Verwendung gelangten, wahrend FLETCHER kugel- 
. f6rmige Teilchen voraussetzte. 

Bei seiner jiingsten Untersuchung der Eiskeimbildung an verschiedenen Metall- 
oxyden hat SERpPoLay [10] fiir ein Kupferoxyd einen ahnlichen Verlauf von ¢(T) 
erhalten, doch erreicht die Kurve schon bei —20°C Werte von g = 0,5. Fiir diese 
Diskrepanz k6nnte einerseits die unterschiedliche Messmethode verantwortlich 
sein, liegen doch die Partikeln bei SERPOLAY auf einer gekiihlten Unterlage, die den 
Eiskeimbildungsmechanismus eventuell auch beeinflusst; der zweite Grund aber 
ist in der verschiedenen Herstellungsmethode zu suchen, womit SERPOLAY auf 
eigentliche Katalysatoren abzielt, die durchwegs sehr porés sind und damit eine 
ausserordentlich grosse spezifische Oberflache aufweisen. 

Bemerkenswert an unseren Messungen ist, dass Stoffe, die auf Grund ihrer 
Schwellentemperatur als gute Eiskeimbildungskerne anzusprechen sind, auch bei 
relativ tiefen Temperaturen eine sehr kleine Aktivitat von nur einigen Promille 
oder Prozent aufweisen. Dies liesse sich damit erklaren, dass die Dichte der fiir 
die Eiskeimbildung giinstigen Stellen auf der Oberflache der Impfsubstanzen recht 
gering ist. Um eine befriedigende theoretische Interpretation des Verlaufs von q(T) 
finden zu kénnen, bedarf es jedoch vergangig noch weiterer experimenteller Unter- 
suchungen. 
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Abschliessend sei dem Leiter des Laboratoriums fiir Atmospharenphysik, Herrn 
Prof. Dr. R. SANGER, fiir die Anregung zu dieser Arbeit und sein stetes Interesse an 
deren Fortgang gedankt. Des Dankes versichert seien auch Herr M. Baup fiir die 
erwahnte Mitwirkung sowie die Herren Prof. Dr. G. Buscu, Prof. Dr. W. EPPRECHT, 
Prof. Dr. F. Laves, Prof. Dr. B. MarINceK und deren Mitarbeiter, die fiir die 
Herstellung und kristallographischen Untersuchungen der Substanzen ihre Appa- 
raturen zur Verfiigung stellten und mit Ratschlagen zur Seite standen. Die finan- 
ziellen Mittel fiir die Durchfiihrung der Untersuchungen sind in verdankenswerter 
Weise von der Schweizerischen Forschungskommission zugesprochen worden. 
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Summary 


The icing nucleability of Cu,O, CuO, Cu,S, Cu,S;, and CuS has been measured 
as a function of the temperature. After a short description of the way in which the 
substances were prepared, the measuring apparatus is explained. Figure 7 shows the 
cloud chamber with a particulary small temperature gradient and the apparatus 
with which a supercooled cloud is first produced in the chamber, before the seeding 
substance from the disperser (Figure 2) is injected into the cloud by means of the 
dosing instrument. The fall-out of ice crystals is observed visually and their number 
estimated. The concentration of the seeding particles is obtained through a sedi- 
mentation tube (Figure 2). The size distribution of the particles is represented 
in Figure 3, and it should be noted that no diameter of less than 0-5 mw occurs. 
Figures 4 and 5 show the activity in terms of icing nucleability q, which is defined 
as the ratio of observed ice crystals to the number of nuclei injected into the cloud. 
A brief comparison is made of these measurements with the theory of FLETCHER [3] 
concerning the dependence of particle size on the threshold temperature, and 
with the experiments of other authors. The only explanation attempted of the 
differences which have been established here has been for the time being in quali- 
tative terms. The threshold temperatures of the substances in question are 
Eee! Cu,O: —4-2°C, CuO: —13-5°C, CuS: —5-1°C, Cu,S.: —o7 C, Colceam 


(Eingegangen: 21, Dezember 1959.) 
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Uber die Trennung von Gasgemischen in laminarer Stromung 


Von JANEZ StRNAD, Ljubljana, Jugoslawien!) 


Es wird eine nur vom theoretischen Standpunkte aus interessante Méglichkeit 
fiir die Trennung von Komponenten eines Gasgemisches betrachtet. 
Fiir ein laminar (Reynoldsche Zahl Re < 2000) durch ein langes Rohr mit 
dem Durchmesser 2 v, str6mendes Gas ist die Geschwindigkeit mit 
y2 
yee BLO) PL = 
(7) = (0) (1— 7) 


gegeben. 
Als Folge der inneren Reibung wird in der Strémung pro Zeiteinheit und Volu- 


meneinheit [1]?) 
2 2 


dv \2 Y 
p=h (+) = 4 wu v?(0) ve 


kinetische Energie in Warme umgewandelt (u Zahigkeit). Wenn die Warmeleitung 
in der Str6mungsrichtung vernachlassigt wird, folgt aus der Warmeleitungsglei- 
chung fiir den stationaren Fall 
WS ae ea b. p ‘ 
YY @ (’ 7) A af Sa 


A Av 
(A Warmeleitfahigkeit). Daraus ergibt sich fiir den Verlauf der absoluten Temperatur 


ce (: ) 


q 
v% 


wobei die Rohrwandtemperatur TJ(v,) als durchwegs konstant angenommen 
wurde. 

Wenn die Stromung hydrodynamisch (Knudsensche Zahl Kn < 0,1) ist, bringt 
die mit T(r) bedingte Thermodiffusion ein Konzentrationsgefalle zustande. In 
einem bindren Gasgemisch, in dem n(v) der Molenbruch der schwereren Kompo- 
nente ist, ergibt sich [2] im stationaren Zustand die Gleichung 


n(r) - (79) [2 + pe v2(0) (1 aie 


1—n(r) = 1—n(7%) 44 T(r) va 


wobei o der Thermodiffusionsfaktor ist. Da im allgemeinen das zweite Glied in der 
grossen Klammer sehr klein ist, kann der Ausdruck in eine Potenzreihe entwickelt 
werden. Wenn man nur die ersten zwei Glieder beibehalt und noch die Konstante 
c = (x w R)/(42M) und die Machsche Zahl Ma = v(0) [M/x R T(r,)}!? einfiihrt, hat 
man 


nip) = nro) — (ya) (1 = m(r9)] ¢¢ Ma? (1 ~ 7) 


(x Adiabatenexponent, M mittleres Molekulargewicht). In den meisten Gemischen 


1) Physikalisches Institut der Universitat. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 246. 
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(c > 0) wird die schwere Komponente am Rande der Stromung angereichert. Der 
Elementareffekt der Trennung e4 = (7p) [1 — n(0)] {n(0) [1 — n(7o)]}-* betragt 


a = COPING. 


Beim Gemisch zweier Isotope mit den Molekulargewichten M, und M, [3] ist o 
von der Gréssenordnung 0,89 (M, — M,)/(M, + M,) und hat e4 somit die Grossen- 
ordnung 0,1 (IM, — M,) Ma?/(M, + M,). Es ware nicht zweckmassig, diesen Effekt 
praktisch auszuniitzen, da andere Trennverfahren*) weit grossere Trenneffekte auf- 
weisen. 

Einen Trenneffekt in laminarer Stromung hat ToLLERT [4] experimentell nach- 
gewiesen. Seine Experimentalbedingungen unterscheiden sich aber von unseren 
Annahmen. 


Tafel von c [5] 


LITERATURVERZEICHNIS 


[1] N. J. Korscurn, I. A. KrpeL und N.W. Rose, Theoretische Hydromechanth — 
Il. Teil (Akademie-Verlag, Berlin 1954), S. 271. : 


[2] S. CHapMAN and T.G.Cow1z1ne, The Mathematical Theory of Non-Uniform 
Gases (Cambridge 1939), S. 244, Gl. (14, 1, 1). : 


[3] R. C. Jones and W. N. Furry, Rev. Mod. Phys. 78, 151 (1946). 
[4] H. Totrert, Z. Elektrochem. 67, 1224 (1957). 


[5] E. H. Kennarp, Kinetic Theory of Gases (McGraw-Hill, New York 1938), S. 180, 
228 
Summary 


The separation effect in a gas mixture laminary and hydrodynamically flowing 
through a cylindrical pipe with constant wall temperature is calculated. 


(Eingegangen: 30. November 1959). 


Cie * 


3) Bei der Diffusion durch pordse Membranen hat zum Beispiel der Elementareffekt der 
Treznung maximal den Wert (M, — M,)/(M, + M,). 
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Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLucee, Band 52: Astrophysik, 111: Das Sonnensystem (Springer-Verlag, 
Berlin 1959). 601 S., 225 Abb.; DM 148.-. 

Von den 600 Seiten des vorliegenden, die Kérper des Sonnensystems behan- 
delnden Bandes sind 362 der Sonne gewidmet. Der einleitende, 79 Seiten umfas- 
sende Abschnitt von GoLpBERG und Pierce itiber die Photosphare gibt eine 
erschopfende Darstellung des kontinuierlichen und des Linienspektrums der 
Sonne, wobei das Schwergewicht auf den Beobachtungen liegt; theoretische 
Untersuchungen, die schon in einem friiheren Band ihre Darstellung gefunden 
haben, werden nur so weit beigezogen, als es fiir die Interpretation und insbesonde- 
re fiir die Erschliessung des Aufbaues der Photosphare notwendig ist. Das 
Kernstiick des Bandes bildet der von DE JAGER verfasste 283 Seiten starke 
Beitrag tiber die Struktur und Dynamik der Sonnenatmosphire. Das Gebiet ist 
so weit gespannt, dass ein einzelner Forscher nur einen kleinen Teil desselben aus 
eigener Tatigkeit kennen kann, im iibrigen auf das Literaturstudium angewiesen 
ist. Trotzdem zeichnet sich der Artikel durch bemerkenswerte Vollstandigkeit 
und Homogenitat aus. Der Autor hat, auch auf ihm ferner liegenden Gebieten, 
mit sicherem Instinkt stets die wesentlichen Punkte hervorgehoben und auch rus- 
sische Arbeiten nicht iibersehen. Wenn man auch dafiir Verstandnis haben kann, 
dass der Autor im wesentlichen nur die Literatur ab 1953 zitiert, so hatte man 
doch die diesbeziigliche Bemerkung lieber am Anfang gesehen, denn es gibt auch 
heute noch Leute, welche einen Artikel vorne zu lesen beginnen; auch hatte da- 
durch manchem Missverstandnis, wo der Leser statt des Namens eines Pioniers 
denjenigen seiner Epigonen findet, vorgebeugt werden k6nnen. Abbildung 47 ist 
voll innerer Widerspriiche, und in der ungliicklichen Abbildung 92 ist durch Um- 
kopieren nach Belieben Schwarz zu Weiss und Weiss zu Schwarz gemacht. 
Trotzdem: alles in allem ein ausgezeichneter Beitrag, der durch seine einfache 
Sprache auch den Lesern nichtenglischer Zunge leicht zuganglich ist. Die nach- 
folgenden kurzen Kapitel behandeln die Atmospharen der Planeten (H.C. 
Urey), den inneren Bau der Planeten (W. C. DE Marcus), die Radioechos vom 
Mond (F. J. Kerr), die Kometen (K. Wurm) und die Meteore (F. L. WHIPPLE). 

Bedauerlich ist, dass der Mond aus dem schén ausgestatteten Band verbannt 
ist, erfreulich, dass derselbe mit nur wenigen Monaten Verspatung auf die neuesten 
darin zitierten Arbeiten erscheinen konnte. M. WALDMEIER 


Atlas of Gamma-Ray Spectra from Radiative Capture of Thermal 
- Neutrons. Von L. D. GrosHeEv, N. Lutsenxo, A.M. DEmMopov und V. I. PELEKHOV. 
Aus dem Russischen iibersetzt von J. B. SyKEs (Pergamon Press, London 1959). 
198 S., 142 Abb.; $20.-. 

Die russischen Autoren haben keine Miihe gescheut, sdmtliche bekannten 
Gammaspektren, die durch Einfang thermischer Neutronen entstanden sind, nach 
Elementen zu ordnen. Fiir dieses Vorhaben sind sie besonders ausgezeichnet, da 
gerade aus dieser Forschungsgruppe viele wertvolle Beitrage herrthren, die aus 
Messungen mit einem Compton-Spektrometer stammen. Das Einfangs-Gammastrah- 
lenspektrum ist in der Regel sehr komplex. Die magischen Kerne mit abgeschlos- 
senen Schalen und einige leichte Kerne bilden eine Ausnahine. Die Kenntnis der 
_ Spektren und besonders der Intensitatsverhaltnisse sind fir viele Arbeiten an 
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einem Reaktor ausschlaggebend. Dieser Atlas ist daher an jedem Reaktor ein un- 
entbehrlicher Ratgeber. Der Kernphysiker findet darin viele Zerfallsschemas, die 
er sonst miihselig in den verschiedensten Handbiichern heraussuchen miisste. Im 
Anhang sind einige Spektren von Konversionselektronen aufgefiihrt, die ebenfalls 
nach Einfang thermischer Neutronen entstehen. Eine sorgfaltige Zusammenstellung 
der Literatur auf diesem Arbeitsgebiet bietet auch dem Spezialisten viel Anregung. 
Als wertvolles Nachschlagewerk diirfte dieser Atlas seine Aufgabe voll und ganz 
erfiillen. P. “STOLE 


Handbook of Probability and Statistics with Tables. Von R.S. Bu- 
RINGTON und D.C. May, Jr. (McGraw-Hill Book Company, New York 
1959). 332S., 53 Fig.; $6.-. 

Dieses schon 1953 durch die Handbook Publishers (Sandusky, Ohio) herausge- 
gebene kleine Handbuch wurde kiirzlich von der McGraw-Hill Book Company ~ 
iibernommen. Da es sich auf die Beschreibung der gebrauchlichsten Methoden der 
mathematischen Statistik beschrankt, hat es seit seinem ersten Erscheinen kaum 
an Bedeutung verloren. Die Verfasser wenden sich sowohl an den Studierenden ~ 
der mathematischen Statistik als auch an den praktisch tatigen Statistiker. Indem 
sie uns einen Uberblick iiber die wichtigsten Zusammenhange geben, berichten — 
sie zugleich alles, was man zur korrekten Anwendung der besprochenen Methoden 
wissen muss. Theorie und Anwendung sind dabei geschickt verkniipft, und die 
mathematischen Begriffe und Beziehungen sind elegant erklart. Sehr wertvoll © 
ware ein ausfiihrlicher Literaturnachweis gewesen. 

Einige Kapiteliiberschriften mégen den Inhalt andeuten: Elementare Wahr- 
scheinlichkeitsrechnung, Verteilungsfunktionen, Priifverfahren und Vertrauens-— 
grenzen, Korrelationsvechnung, Zeitserien, Streuungszerlegung, Fortschreitende Pro- 
benahme, Qualitatskontrolle. Neben den tiblichen Tafeln fiir +, F- und y?-Test 
finden wir eine Anzahl niitzlicher Tabellen. Spezielle Erwahnung verdienen die — 
ausfihrlichen Tafeln fiir Binomial-, Poisson- und Normalverteilung. 4 
P.SCHMID © 


Basic Physics of Atoms and Molecules. Von U. Fano und Mrs. L. Fano 
(John Wiley & Sons, New York 1959). 414 S., 86 Fig.; $10.-. , 

Das vorliegende Buch postuliert in der Einleitung eine Liicke im vorhandenen — 
Schrifttum, die es nun ausfiillen will. Die Liicke besteht im Mangel an Werken, i 
welche eine Einfiihrung und ein Versténdnis fiir die Mikrophysik und ihren Zusam- 
menhang mit den gewohnten Eigenschaften materoskopischer Objekte geben. 

_ Das Haupthindernis fiir das Eindringen in das Gebiet der Quantentheorie sehen — 
die Autoren wohl mit Recht in dem dazu erforderlichen mathematischen Apparat. 

Unter diesen Voraussetzungen gehen die Autoren zum Teil eigenwillige Wege. 
So werden wichtige quantenmechanische Tatsachen zundchst fiir diskrete Situa-_ 
tionen (Quantisierung des Drehmomentes) illustriert. Nun ist das Diskrete allerdings } 
begrifflich einfacher als das Kontinuierliche, doch folgt daraus nicht notwendig, — 
dass es zur Einfithrung besser geeignet ist. 

Die Autoren machen sich ihre Aufgabe nicht leicht, indem sie ihre eigenen Wege : 
gehen. Das Buch ist durchaus sorgfaltig geschrieben und behandelt neben dem 
itblichen Stoff der Atomphysik zur Hauptsache noch die chemische Bindung. Ein : 
Reihe von Anhangen sind der Erganzung der mathematischen Hilfsmittel gewidmet. 

Das Werk richtet sich in erster Linie an Fachleute aus den Nachbarwissenschaf- 
ten zur Physik, ausgenommen an Mathematiker, also etwa an Chemiker und 


Biclogen. R Jost 


Voranzeige 


Mathematisthe Reithe, Band 25 und 26 
Sammlung «Lehrbiicher und Monographien aus dem Gebiete der exakten Wissenschaften» 


ACZEL 


‘Vorlesungen iiber Funktionalgleichungen 
und ihre Anwendungen 


Von Prof. Dr. J. Aczél, Professor an der Universitat Debrecen (Ungarn). Ca. 320 Seiten mit 
14 Figuren. Préis gebunden ca. Fr./DM 40.— 

Inhalt: Inhaltsverzeichnis. Einleitung. — Erster Teil: Gleichungen fiir Funktionen von einer 
Veranderlichen: 1. Gleichungen, in denen Veranderliche nicht nur unter Funktionenzeichen _ 
vorkommen. 2. Gleichungen mit lauterinneren Veranderlichen. 3. Gleichungen mit mehreren 
unbekannten Funktionen. 4. Zuriickfiihrung. auf Differential- und Integralgleichungen. — 
Zwetter Tetl: Gleichungen fiir Funktionen mehrerer Veranderlicher: 5. Einfache Gleichungen. 
6. Zusammengesetzte Gleichungen. 7. Gleichungen mit mehreren unbekannten Funktionen 
mehrerer Veranderlicher. Zuriickfiihrung auf partielle Differentialgleichungen. 8. Vektoren- 
und Matrizengleichungen. — Schlussbemerkungen. Einige ungeléste Probleme. — Literatur- 
; verzeichnis. 

Das vorliegende Buch ist das erste, das ein grosses Gebiet der Theorie der Funktionalglei- 
chungen zusammenfassend bearbeitet. Die grossen Anwendungsgebiete, insbesondere Wahr- 
scheinlichkeitstheorie (nichteuklidische) Geometrie und Mechanik, werden auch eingehend 
behandelt. Das Literaturverzeichnis, das mehr als 1000 Arbeiten aufzahlt, will gewisser- 

massen auch als Bibliographie dieses Gebietes dienen. 


KOWALSKY 
Topologische Raume 


Von Prof. Dr. H.-J. Kowatsxy, Professor an der Universitat Erlangen 
Ca. 240 Seiten, Preis ca. Fr./DM 35.— 


Inhali: 1. Grundlagen: Mengen, Verbande, Filter. Il. Topologische Raume: Zusammenhang 
zwischen Topologie und Grenzwert, Grundbegriffe, Trennungseigenschaften, Machtigkeits- 
bedingungen. III. Kompaktheit und Zusammenhang (einschliesslich lokaler Begriffe und 
parakompakter Raume). IV. Abbildungen: Stetige, offene und abgeschlossene Abbildungen, 
Homéomorphie, vollstandig regulare Raume, Quotienten-, Produkt-, Summen- und Ab- 
bildungsraume. V. Erweiterung und Kennzeichnung topologischer Raume: Allgemeines 
Erweiterungsprinzip, Kompaktifizierung, Einbettungs- und Darstellungssatze. VI. Metrische 
und uniforme Raume: Metrisierung, Gleichmassigkeit, Vervollstandigung. VII. Topologi- 
sche Gruppen, Anwendungen: Grundbegriffe aus der Theorie der topologischen Gruppen, 
Ringe, Kérper und Vektorraume, Approximationssatz von Stone-Weierstrass, induktiver 
und inverser Limes. 


Das vorliegende Lehrbuch soll in die mengentheoretische Topologie einfiihren. Ausgangs- , 
punkt ist eine Verallgemeinerung des Grenzwertbegriffs, die einen natiizlichen Zugang zu 
den topologischen Raumen vermittelt. Kennzeichnend fir die Darstellung ist die weit- 
’ gehende Verwendung der Filter, die nicht nur eine tibersichtliche Schreibweise und Begriffs- 
bildung, sondern aiich eine formale Beweisfiithrung erméglichen. Die Stoffauswahl wird 
durch Hinweise und Aufgaben am Ende eines jeden Paragraphen erganzt. 
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Vorlesungen iiber Differential- und Integralrechnung | 
In drei Banden, von A. OSTROWSKI, Professor an der Universitat Basel. 


Zum Gebrauch bei akademischen Vortragen sowie zum Selbststudium 


Mathematische Reihe, Band 4, 5 und 7 — Sammlung «Lehrbucher und Monographien aus _ ie 
dem Gebiete der exakten Wissenschaften» i 


Soeben erschienen: 
Band I: Funktionen einer Variablen 
Zweite, neubearbeitete Auflage. 330 Seiten mit 47 Figuren. Preis Fr./DM 35.— (1960) 


In der neuen Auflage haben die beiden ersten Bande gewisse strukturelle Anderungen er- | 
fahren. Um den ersten Band, mit dem erfahrungsgemass fiir manche Studierende der Lehr- 

"gang der Infinitesimalrechnung abgeschlossen wird, in sich abgerundet zu machen, wurden 

aus dem zweiten Band die Elemente der Reihenlehre und einige Grundtatsachen der 
Differentialgeometrie heriibergenommen. ieee as 

Inhalt; Wesen der Mathematik. I. Grundbegriffe: § 1. Die Grundeigenschaften reeller Zahlen. ! 

§ 2. Erste Folgerungen aus den Grundeigenschaften reeller Zahlen. § 3. Weitere Folgerungen 

aus den Grundeigenschaften. § 4. Der Funktionsbegriff. II. Grenzwerte: § 5. Nullfolgen. 

§ 6. Grenzwerte von Zahlenfolgen. § 7. Spezielle Satze tiber Konvergenz von Zahlenfolgen. 

§ 8. Unendliche Reihen. § 9. Grenzwerte von Funktionen eines stetigen Arguments. JI. 
Stetige Funktionen einer Variablen und bestimmte Integrale: § 10. Stetige Funktionen. § 11. 
Trigonometrische Funktionen. § 12. Definition des bestimmten Integrals. § 13. Elementare 

* Eigenschaften des bestimmten Integrals. IV. Der Begriff der Ableitung und die Fundamen- 
talsdtze der Infinitesimalrechnung: § 14. Die Ableitung. § 15. Der Zusammenhang zwischen. 

der Ableitung und dem Differenzenquotienten einer Funktion. § 16. Die Fundamentalsatze 

der Infinitesimalrechnung. V. Die Technik des Differenzterens: § 17. Ableitungen rationaler 
Verbindungen gegebener Funktionen. § 18. Umkehrung monotoner Funktionen. § 19, Die 
Kettenregel und ihre Anwendungen. VI. Die Technik des Integrierens: § 20. Integralum- 
formung durch partielle Integration und Variablensubstitution. § 21. Logarithmus und 
Exponentialfunktion. § 22. Anwendungen der Logarithmus- und der Exponentialfunktion. | 

§ 23. Integration rationaler Funktionen. VII. Erste Anwendungen der Differentialrechnung | 

auf die Geometrie und die Funktionen-Diskussion: § 24. Anwendungen der ersten Ableitung. 

; Monotonie. § 25. Anwendungen der ersten Ableitung. Extrema. Unbestimmte Ausdriicke. 
ieee yt § 26. Die zweite Ableitung. Der Funktionsverlauf im Grossen. § 27. Darstellung von Kurven | 
in der Ebene und im Raum. § 28. Bogenlange. § 29. Entwicklung des Logarithmus und des | _ 
Arcustangens. § 30. Die Taylorsche Formel. § 31. Entwicklungen von e”, sinx und cos%. 
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Band II: Differentialrechnung auf dem Gebiete mehrerer Variablen 
484 Seiten mit 55 Figuren. Leinen Fr./DM 67.- (1951) 
Inhalt: Unendliche Mengen — Funktionen auf Mengen — Unendliche Folgen und Reihen — 
Erganzungen zur Differentialrechnung — Anwendungen der Differentialrechnung auf die 
Analysis — Numerische Rechenmethoden — Bogenlange — Ebene Kurven — Raumkurven 
"und Flachen. q 


Band III: Integralrechnung auf dem Gebiete mehrerer Variablen | 

475 Seiten mit 86 Figuren. Ganzleinen Fr./DM 75.— (1954) i 

Inhalt: Erganzungen zur Integrationstechnik — Begriff des mehrfachen Integrals - Berech- i 
nung mehrfacher Integrale - Anwendungen mehrfacher Integrale — Einfache uneigentliche | 
Integrale — Mehrfache uneigentliche Integrale und die Gammafunktion — Fourier-Reihe | 
und Fourier-Integrale. iy 
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